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\  Abstract 

-V 

Tha  purpose  of  this  study  was  to  investigate  the 
perforeance  characteristics  of  a  synchronization  systea  that 
uses  an  auxiliary  matched  filter  channel  to  speed  accpiisition 
tine.  The  study  entailed  a  quantitative  analysis  and 
simulation  of  the  matched  filter  channel.  Following  a 
presentation  on  computer  simulation  of  communication  systems, 
the  investigation  took  the  following  approach:  (1)  Perform  a 
quantitative  analysis  to  characterize  mean  acquisition  time 
under  varying  conditions  of  matched  filter  length,  noise 
levels,  and  number  of  users  in  the  system.  (2)  Build  and 
validate  a  simulation  model  of  the  synchronization  system 
using  the  Block  Oriented  Systems  Simulator  (BOSS) .  (3)  Test 

the  matched  filter  channel  against  the  analytic  expectations 
for  false  alarm  performance  via  the  simulation  model. 

The  quantitative  results  indicate  that  the  matched 
filter  channel  improves  acquisition  time.  For  signal-to> 
noise  ratios  above  OdB  and  a  particular  matched  filter 
length,  the  mean  acquisition  time  was  directly  a  function  of 
the  number  of  users  in  the  system.  For  a  matched  filter  of 
length  of  128  stages  and  64  users  in\^e  system,  the  matched 
filter  channel  provides  a  gain  in  acquisition  time  of  156 
over  the  active  correlator  device  operating  alone.  As 
signal-to-noise  ratios  decreased  below  OdB,  the  variation  in 


■ean  acquisition  tise  caused  by  the  number  of  user 3  declines 
and  the  additive  white  Gaussian  noise  becomes  the  dominant 
influence.  Mean  acquisition  time  asymptotically  approaches  a 
doubling  in  time  for  each  1.5  dB  decrease  in  the  signal-to~ 
noise  ratio. 

The  simula-cion  results  confirmed  the  analytic 
expectations  for  false  alarm  performance.  The  simulated 
probability  of  false  alarm  (Pj^)  was  restricted  to  lo”^.  All 
single-pulse  detection  requirements  at  the  specified  were 
favorably  confirmed.  Multiple-pulse  detection  requirements 
were  favorably  validated  with  shorter  simulations. 

Exceptions  to  the  multiple-pulse  detection  results  were  noted 
due  to  the  high  correlation  noise  produced  by  the  selected 
synchronization  codes. 

The  BOSS  provided  a  highly  flexible  simulation 
environment.  The  system  model  included  flexibilities  to 
adjust  matched  filter  lengths,  noise  levels,  and  threshold 
levels.  Confirmation  of  the  analytic  expectations  provides  a 
proven  simulation  model  to  accurately  simulate  dynamic 
applications  of  system. 


ANALYSIS  AND  SIMULATION 
OF  A 

PSEUDONOISE  SYNCHRONIZATION  SYSTEM 

Introduction 

Computer  simulation  provides  a  flexible  and  accurate 
method  for  performance  evaluation  and  trade-off  analysis  in 
the  design  of  digital  communication  systems.  Recent  advances 
in  simulation  techniques  and  advances  in  computer  software 
and  hardware  technologies  have  merged  to  provide  intelligent, 
user  friendly  and  flexible  simulation  environments. 

Recently,  there  has  been  a  growing  interest  in  the  use  of 
simulation-based  techniques  for  the  analysis  and  design  of 
communication  systems  (1:36.1.1). 

In  order  to  design  and  analyze  complex  spread  spectrum 
systems  effectively  and  continue  to  meet  the  stringent 
performance  req[uirements ,  design  techniques  have  evolved  that 
combine  traditional  analysis  with  simulation.  In  fact, 
simulation  has  almost  become  the  standard  tool  for 
communication  system  analysis  and  design  (2:1). 

Problem  Statement 

In  spread  spectrum  communication  systems,  a  basic  task 
to  be  performed  at  the  receiving  end  is  the  synchronization 
of  the  pseudonoise  (PN)  signal  generated  locally  at  the 
receiver  with  the  PN  signal  contained  in  the  receiver  input 
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signal.  The  problem  addressed  in  this  thesis  is  to  analyze 
and  simulate  a  synchronization  scheme  that  uses  an  auxiliary 
matched  filter  channel  to  speed  synchronization  and  evaluate 
the  initial  acquisition  characteristics.  Such  a 
synchronization  system  has  been  described  by  Smirnov  (3) . 

Scope 

The  Block-Oriented  Systems  Simulator  (BOSS)  will  be  used 
to  simulate  the  synchronization  system.  The  simulation 
models  will  include  flexibilities  to  vary  the  matched  filter 
length,  adjust  threshold  detection  levels,  and  specify  the 
number  and  noise  levels  of  transmitted  signals.  Once  the 
system  model  is  validated,  further  simulations  will  exercise 
the  system  parameters  to  fully  characterize  acquisition 
performance  and  confirm  analytic  expectations.  The  key 
performance  factor  is  the  mean  acquisition  time  as  impacted 
by  the  matched  filter  length  and  the  number  of  users  in  a 
code  division  multiple-access  (CDMA)  environment. 

ttatgrialg 

•  BOSS  User's  Manual  (4) 

■  BOSS  Software  Version:  ST*AR  1.2 

•  DEC  VAXstation  II 

Approach 

The  overall  design  philosophy  is  to  implement  the 
synchronization  system  incrementally  in  the  simulation 
environment.  Each  system  module  is  implemented,  tested 
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against  analytic  expectations,  and  stored  in  the  simulation 
module  library.  The  complete  synchronization  system  is  then 
constructed  and  analyzed  with  a  single  transmitted  signal. 
Subsequent  simulation  and  analysis  is  then  performed  for  a 
multiple-access  environment. 

Thggjg  Qrqanizfltigp 

In  the  next  chapter,  computer  simulation  is  discussed. 

In  Chapter  III,  the  system  under  simulation  is  presented,  the 
synchronization  method  is  introduced,  and  the  necessary 
assumptions  are  justified.  An  analysis  concerning  the 
synchronization  time  and  the  impact  of  the  matched  filter 
length  and  the  number  of  users  is  given,  chapter  IV  is 
devoted  to  the  modeling  of  the  necessary  system  modules  using 
the  BOSS.  The  chapter  is  concluded  with  the  integration  of 
the  modules  for  the  complete  system  representation.  In 
Chapter  V,  the  objective  and  method  of  simulation  is  detailed 
and  the  simulation  results  are  presented.  Finally,  in 
Chapter  VI,  conclusions  on  the  analysis  and  the  effectiveness 
of  the  simulation  are  given  with  recommendations  for  further 


research . 


Computer  simulation  has  almost  become  the  standard  tool 
for  communication  system  analysis  and  design  (2:1). 

Simulation  offers  flexibility  and  insight  into  system 
performance.  Although  computer  simulation  has  distinct 
advantages,  the  approach  tends  to  be  distrusted  by  the 
proponents  of  analytical  modeling  (5:89).  Simulation  yields 
answers  with  statistical  variations,  not  hard  results. 
Analytic  techniques,  however,  are  not  adequate  to  predict  the 
performance  of  systems  operating  in  complex  environments. 

The  presence  of  nonlinearities,  noise,  and  other  degradations 
make  the  performance  evaluation  a  tedious  and  often 
impossible  task  (6:1612).  Techniques  have  emerged  that 
combine  simulation  and  analysis.  "The  traditional  dichotomy 
between  analysis  and  simulation  has  now  been  blurred"  (2:1). 

A  clear  understanding  of  simulation  techniques  and 
methodology  becomes  Indispensable  as  stringent  performance 
requirements  drive  engineers  to  effectively  design  and 
analyze  increasingly  complex  systems. 

This  chapter  is  a  short  primer  on  computer  simulation 
techniques  in  regards  to  communication  system  simulation. 

The  first  section  will  look  at  the  basic  operations  in 
communication  system  modeling  and  introduce  simulation  as  a 
technique  to  estimate  system  performance.  The  next  section 
will  extend  quantitative  performance  evaluation  via  the  Monte 
carlo  method  as  the  primary  technique  of  interest.  The 


critical  features  of  a  comaunication  system  simulation 
package  will  be  examined  and  finally,  the  methodology  of  the 
Block  Oriented  System  Simulator  (BOSS)  will  be  explored. 

SlBUlation  Approach 

There  are  four  steps  fundamental  to  the  simulation  of  an 
elementary  communication  channel:  1)  create  and  store  a 
sequence  of  samples  representing  the  modulated  signal; 

2)  filter  the  sequence  of  samples;  3)  modify  the  filtered 
signal  according  to  channel  nonlinearities;  and  4)  demodulate 
the  data  record  to  observe  the  effect  of  the  channel  (5:89). 
These  four  steps  will  be  examined  in  terms  of  the  elementary 
channel  shown  in  Figure  1. 


Data 

Data 

Modulator 

- 

Linear 

Filter 

• 

Non- 

Linearity 

Demodulator 

- ► 

- ► 

Figure  l.  Elementary  Communication  Channel  (5:90) 


Signal  Generation.  The  first  step  is  to  create  a  set  of 
sampled  signals.  Modulated  signals  are  generated  as  sampled 
versions  of  the  real  signal.  For  a  sampling  frequency 
f  >  1/dt  the  complex  envelope  representation  of  the 

S 

k-th  sample  is 


s^(k  dt)  -  A(k  dt)yr  exp[j2s(f^f^)k+j©(k  dt)]  (1) 
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where  A(k'dt)  represents  the  amplitude  modulation,  and 
6(k-dt)  represents  the  phase  modulation  of  the  signal,  with 
power  P  and  relative  carrier  frequency  (5:90). 

The  prime  consideration  in  selecting  f^  is  the  total 
bandwidth  about  the  channel  carrier  frequency.  In  spread 
spectrum  systems,  this  bandwidth  is  determined  by  the 
frequency  of  the  spreading  waveform.  The  sampling  rate  is 
expressed  as  a  multiple  N  of  the  spreading  code  chip  rate: 


N/Tc 

(2) 

VN 

(3) 

The  number  of  samples  per  chip  period  is  usually  selected  as 
an  even  integer  between  8  and  16.  For  values  less  than  8 
samples  per  chip  period,  accuracy  is  sacrificed  at  the 
folding  frequencies  due  to  aliasing  (5:90). 

Filtering  of  Signals.  Filtering  tends  to  be  the  most 
CPU-intensive  activity  in  channel  simulations.  There  are  two 
popular  techniques  to  simulate  filtering:  1)  multiply  the 
discrete  Fourier  transform  (OFT)  of  the  signal  by  the  filter 
function  and  return  to  the  time  domain  by  implementing  the 
Inverse  DFT;  and  2)  operate  on  the  sampled  signal  using  a 
digital  filtering  algorithm  (5:91) . 

In  the  first  technique,  the  signal  data  record  is 
segmented  into  conveniently  sized  blocks  (512,  1024,  etc.) 
of  data.  A  DFT  of  these  smaller  blocks  of  data  is  then 
synthesized  and  multiplied  by  the  filter  function.  Once  this 
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operation  is  performed  on  the  total  data  record,  an  inverse 
OFT  is  performed  to  yield  the  filtered  time  domain  waveform. 

In  many  applications,  the  filters  are  classical  types 
such  a  Butterworth,  Chebyshev,  or  Bessel.  To  simulate  these 
filters  and  those  in  which  pole^zero  descriptions  are 
available,  it  is  convenient  to  synthesize  a  recursive  digital 
filter  to  implement  a  time  domain  filtering  algorithm. 

Channel  Nonl Inearit ies .  When  the  channel  model  contains 
nonlinearities,  analytical  approaches  become  complex,  and 
simulation  is  often  the  only  practical  approach  to 
performance  evaluation  (5:92).  The  general  approach  to 
simulating  channel  nonlinearities  involves  the  use  of 
measured  data  of  representative  nonlinear  devices.  The 
measured  data  is  used  to  generate  a  function  for  mapping 
input  to  output  values.  For  abrupt  nonlinearities,  table 
lookup  procedures  are  used  to  define  nonlinear 
characteristics . 

Receiver  Structures.  The  first  element  of  the  receiver 
structure  is  the  receiver  front-end  filter.  At  this  point, 
the  method  of  adding  thermal  noise  at  the  input  to  this 
filter  distinguishes  two  approaches  to  channel  simulation. 

In  the  hybrid  simulation/analysis  approach,  the  simulation 
is  used  only  to  obtain  the  statistics  of  all  other  sources 
and  interference.  The  effect  of  thermal  noise  is  then  added 
analytically  (5:93). 

In  the  simulation  approach  of  interest,  thermal  noise  is 
added  at  appropriate  points  in  the  channel.  This  approach. 
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direct  channel  simulation,  is  the  only  choice  when  the 
demodulation  and  detection  process  in  nonlinear  (5:93).  The 
demodulation  and  detection  process  is  followed  by 
quantitative  evaluation  where  bit  error  rate  (BER)  is  the 
most  often  used  performance  characteristic  for  digital 
communication  systems  (7:153). 

The  demodulation  and  detection  process  output  is 
compared  to  a  decision  threshold.  The  threshold  comparison 
depends  on  the  conditional  probability  density  functions 
p(z|s^)  and  p(z|sq).  Figure  2  shows  representative  densities 
that  in  general,  need  not  be  identical  (7:154).  An  error 
will  occur  when  either:  1)  a  "zero"  is  sent  and  the  channel 
noise  results  in  z(T)  being  greater  than  the  threshold  0^  or 
2)  a  "one"  is  sent  and  the  channel  noise  results  in  z(T) 
being  less  than  0^.  The  probabilities  of  these  errors  are 
given  by 


Prob{error I s^ ) 


Pi 


f.0O 

p(z|sj^)dz 

-« 


and 


Prob(error I Sq) 


Pq 


r«o 

p(z|sQ)dz 

0o 


(4) 


(5) 


The  total  probability  of  error  is  P{s^)  p^  P{Sq)  Pq  where 
P{s^}  and  P{Sq}  are  the  a  priori  probabilities  of  symbols 
"one"  and  "zero"  respectively. 

Techniques  for  estimating  the  BER  in  the  direct 
simulation  approach  include:  1)  Monte  Carlo  simulation, 
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Figure  2.  Conditional  Probability  Density  Functions 


2)  importance  sampling  (modified  Monte  Carlo  simulation) , 

3)  extreme  value  theory,  and  4)  tail  extrapolation.  The  four 
techniques  differ  in  the  way  the  probability  density 
functions  are  handled.  The  Monte  Carlo  method  is  the  most 
general  and  uses  an  empirical  determination  of  the 
distribution  functions  evaluated  at  a  single  point. 

Importance  sampling  biases  the  distribution  functions  to 
cause  errors  to  occur  more  often.  The  BER  is  then  found  by 
unbiasing  the  results.  This  method  can  substantially  reduce 
the  simulation  time.  The  extreme  value  theory  and  tail 
extrapolation  operate  under  the  principle  that  the 
distribution  of  the  extreme>value  (large  in  magnitude) 
samples  is  the  only  information  necessary  to  calculate  the 
BER  (7:154-155).  Of  these  techniques,  only  the  Monte  Carlo 
method  will  be  discussed.  The  BOSS  in  essence  uses  the  Monte 
Carlo  method. 
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KsntA  Cacifi  siffwlfltion 

Basically,  the  Monte  Carlo  aethod  is  an  error  counting 
technique.  In  other  words,  the  BER  can  be  estimated  by 
observing  the  number  of  error  occurrences.  Suppose  that  a 
"zero"  is  sent.  Then  the  conditional  probability  of  error. 
Eg  (5),  can  be  expressed  as  (7:156): 


Po 


.00 

feo(*)P(2|«o)d* 


(6) 


where  h^Cz)  is  basically  an  error  detector: 


,  z  i  0^ 

0,  Z  <  0^ 


(7) 


Recognizing  that  Eq  (6)  expresses  the  expected  value  of 
h^Cz),  the  conditional  probability  of  error  can  be  estimated 
with  a  sample  mean  of  h^Cz).  In  other  words,  p^  >  ECh^Cz)] 
can  be  estimated  by  counting  the  number  of  error  occurences 
over  N  observations  as  given  by: 


^0 


(8) 


The  above  equation  defines  the  Monte  Carlo  method.  An 
estimate  of  the  BER  is  given  by  counting  the  error  occurences 
over  a  set  of  observations.  Simplifying  Eq  (8)  for  "n”  bit 
errors  out  of  "N"  bits  processed  by  the  system  gives  an 
unbiased  Monte  Carlo  estimate  of  the  BER: 
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n 


(9) 


P  •  - 

N 

In  th*  Halt,  as  N  approaches  infinity,  "j^"  will  converge  to 
the  true  value  "p"  (7:156).  Obviously,  it  would  take  an 
infinite  amount  of  computer  time  for  the  estimate  to  converge 
to  the  true  value.  However,  for  a  finite  number  of 
observations,  the  reliability  of  the  estimate  can  be 
quantified  in  terms  of  confidence  levels.  The  confidence 
level  is  often  specified  as 

Problhj  <  p  <  hj^)  »  (1-a)  (10) 

where  (l-a)  is  the  probability  that  the  true  value  is 
bracketed  by  h^  and  h^. 

The  two'sided  confidence  interval  is  plotted  in  Figure  3 

for  confidence  levels  of  90  percent,  95  percent,  and  99 

percent.  The  "rule  of  thumb"  is  that  the  number  of 

observations  should  be  on  the  order  of  10/P_  where  P_ 

e  e 

represents  the  probability  of  error.  This  coincides  with  a 

k+1 

vertical  slice  at  N  »  10  on  Figure  3.  For  a  mildly 
distorted  communication  system,  this  level  of  uncertainty 
corresponds  to  ±0.5  dB  on  the  signal -to-noise  ratio  (Ej^/N^) 
scale  of  a  typical  BER  curve.  This  level  of  uncertainty  is 
generally  considered  reasonable  (7:157). 

When  simulating  a  system,  a  number  of  BER  estimates  are 
made  at  discretely  spaced  values  of  A  best  fit  curve 

is  taken  to  represent  the  characteristic  behavior  of  the 
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Figure  3.  BER  Confidence  Bands  (7:158) 


system.  This  curve  will  take  into  account  the  relative 
reliability  of  the  individual  points  and  provide  a  somewhat 
better  confidence  than  predicted  from  the  single  point 
confidence  Intervals  of  Figure  3. 

Software  Structure 

A  number  of  software  packages  have  been  developed  for 
simulation-based  analysis  and  design  of  communication 
systems.  Most  of  these  software  packages  were  developed 
during  the  seventies  and  their  design  and  implementations  are 
based  on  older  technologies.  Advances  in  computer  hardware 
technology,  software  technology  and  simulation  techniques 
offer  a  totally  new  computing  environment  (4:36.1.1). 


The  main  goal  in  developing  a  simulation  software 
package  is  to  provide  an  integrated  computing  environment 
that  will  allow  a  communication  system  engineer  to  design 
models,  configure  and  execute  simulations,  review  the  results 
and  perform  design  iterations.  While  there  is  no  industry 
standard,  a  software  simulation  package  should  include  a 
modular  structure  and  a  user-friendly  interface  (8:2). 

A  modular  structure  provides  maximum  flexibility  in  the 
design  and  analysis  of  communication  systems.  Most  of  the 
software  packages  now  in  use  are  made  up  of  four  major 
components:  a  system  configurator,  a  model  library,  a 
simulation  exerciser,  and  a  post-processor  (9:324). 

The  system  configurator  is  used  to  select  a  set  of 
pre-programmed  models  of  functional  blocks  from  the  model 
library.  These  functional  blocks  are  then  connected  in  the 
desired  topology  as  specified  by  the  block  diagram  of  the 
system  being  designed.  Parameters  of  various  functional 
blocks  are  specified  either  when  the  system  is  being 
configured  or  during  the  executi)on  of  the  simulation.  The 
simulation  exercisor  supervises  the  execution  of  the 
simulation  (9:324). 

Time  histories  of  events  or  waveforms  at  various  points 
in  the  system  are  generated  and  stored  during  the  simulation. 
These  time  histories  or  waveforms  are  examined  by  the 
post-processor  at  the  conclusion  of  the  simulation  run  and 
performance  measures  are  computed  from  the  time  histories. 

The  post-processor  also  provides  graphic  displays  of 
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wavefons,  spectra,  and  histograns  of  signals  at  various 
points  in  the  system.  At  the  end  of  a  simulation,  the  design 
engineer  uses  the  post^processor  output  to  verify  if  the 
performance  requirements  and  design  constraints  are  met.  If 
design  objectives  are  not  met,  several  iterations  are  made 
until  a  suitable  design  is  found  (9:324). 

Several  software  pac)cages  are  presently  available  for 
simulating  communication  systems.  The  state-of-the-art 
features  and  capabilities  of  the  Block-Oriented  Systems 
Simulator  (BOSS)  are  ideally  suited  for  simulating  spread 
spectrum  systems. 

BOSS  Methodology 

The  Block-Oriented  Systems  Simulator  (BOSS)  is  an 
interactive  environment  for  time-domain  (waveform  level) 
analysis  and  design  of  communication  systems.  The  user 
interface  controls  all  design  capabilities  through  the  use  of 
interactive  graphics.  Capabilities  are  provided  to  design 
models  using  block  diagrams,  configure  and  execute  a 
simulation,  review  the  results  of  the  simulation,  and  perform 
design  iterations  (1:36.1.1,36.1.6). 

The  BOSS  provides  for  simulation  of  any  system  which  can 
be  represented  in  block  diagram  form.  Complex  communication 
systems  are  synthesized  in  a  hierarchical  fashion  using 
simple  building  blocks  at  the  lowest  levels,  and  increasingly 
complex  modules  and  subsystems  at  subsequent  higher  levels. 
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Each  block  of  the  simulation  model  performs  a  signal 
processing  operation  (4:11,13). 

During  the  initial  model  construction  phase,  BOSS 
requires  the  designer  to  provide  complete  descriptions  of 
signals  and  parameters  such  as  data  types,  ranges,  and 
default  values.  Error  checking  is  automatically  performed  by 
BOSS  on  all  input  data.  The  use  of  menus  for  entering  data, 
type  and  range  checking,  and  other  consistency  checks  reduce 
input  errors  and  minimize  confusion  (4:47). 

Displays  of  simulation  results  are  provided  by  the  Post 
Processor.  Multiple  windows  for  display  as  well  as  multiple 
displays  within  a  window  are  standard  features.  Automatic 
scaling  of  plots,  labeling,  pan  and  zoom  are  all  built-in 
features  (4:13). 

BOSS  allows  flexible  interactive  simulation  of  a  variety 
of  design  problems.  This  flexibility  is  fundamental  for 
performance  evaluation  and  trade-off  analysis  in  the 
simulation  and  design  of  digital  communication  systems. 


ill.  Pr9>?lgB  Analysis 


The  purpose  of  this  chapter  is  to  present  in  detail  the 
systen  to  be  simulated.  Prefaced  with  a  background  on  the 
synchronization  of  spread  spectrum  systems,  an  outline  of  the 
general  concepts  of  signal  flow  and  operation  of  the  circuit 
of  interest  is  presented.  An  analytic  evaluation  of  the 
system  is  then  performed  to  present  the  performance 
characteristics  and  to  facilitate  the  proper  structure  and 
evaluation  of  the  simulation  effort. 

Synghr<?ni£a.ti<?n  ol  Spcead  spectrum  gys-t^eas 

For  military  radio  communications,  communication 
security  consists  of  preventing  signals  from  being  detected 
or  jammed  by  an  adversary.  Spread  spectrum  signals  meet  both 
of  these  goals  and  also  offer  the  possibility  of  multiple'- 
access  and  range  measurement.  To  qualify  as  a  spread 
spectrum  system,  a  communications  link  must  meet  two 
criteria.  First,  the  link  must  transmit  signals  using  a  much 
wider  bandwidth  than  is  required  to  send  the  information  of 
interest.  Second,  an  independent  function  other  than  the 
information  being  sent  is  used  to  determine  the  resulting 
modulated  bandwidth.  The  information  independent  technique 
used  to  spread  the  bandwidth  of  a  transmitted  signal  is  also 
used  to  despread  the  received  signal  (10:855). 

The  two  primary  techniques  used  to  accomplish  the  signal 
bandwidth  spreading  and  subsequent  despreading  are  direct 
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sequence  sodulation  and  frequency  hopping.  Direct  sequence 
aodulation  causes  rapid  phase  transitions  in  the  carrier 
frequency  to  spread  the  overall  signal.  The  individual 
transitions  are  randomly  controlled  by  a  pseudorandomly 
generated  sequence.  When  frequency  hopping  is  employed,  a 
pseudorandom  or  pseudonoise  (PN)  code  is  also  used,  but  in 
this  case  the  code  is  used  to  shift  the  entire  carrier 
frequency  and  not  just  its  phase.  In  both  cases,  the 
objective  is  to  randomize  the  modulated  wavefoirm  so  the 
transmitted  bandwidth  depends  almost  entirely  on  the 
pseudorandom  code  used.  If  properly  coded,  the  transmitted 
signal  is  very  difficult  to  acquire  and  despread  without 
knowledge  of  the  code. 

In  order  to  despread  the  received  signal,  the  spread 
spectrum  receiver  generates  a  replica  of  the  PN  code  and  then 
synchronizes  this  local  PN  sequence  to  the  one  used  to  spread 
the  incoming  signal.  Multiplication  of  the  incoming  signal 
by  the  synchronized  local  PN  sequence  produces  the  desired 
despreading  process.  The  problem  of  synchronization  is  that 
of  aligning  the  receiver's  locally  generated  waveform  with 
the  spreading  waveform  superimposed  on  the  incoming  signal 
(11:3). 

Classically,  the  synchronization  process  is  accomplished 
in  two  phases.  Initially,  the  PN  code  acquisition  phase 
brings  the  two  codes  into  coarse  alignment.  Once  the 
incoming  PN  code  has  been  acquired,  the  PN  tracking  phase 
maintains  waveform  alignment  with  a  closed  loop  operation. 


There  are  two  widely  used  techniqpies  for  detection  of 
known  input  signals.  The  first  is  based  on  active 
correlation  where  for  acquisition,  the  local  PN  sequence 
generator  is  clocked  in  such  a  way  that  its  sequence  is 
shifted  continuously  or  in  discrete  steps  relative  to  the 
sequence  contained  in  the  received  signal.  The  instants  of 
synchronism  are  detected  by  the  acquisition  circuit  which  in 
turn  activates  the  control  loop  provided  for  tracking.  The 
second  technique  is  based  on  the  principle  of  passive 
correlation.  The  acquisition  circuit  has  the  task  of 
detecting  the  instant  at  which  the  received  PN  sequence 
attains  a  predetermined  state.  In  the  case  of  passive 
correlation,  the  basic  component  of  the  acquisition  circuit 
is  typically  a  matched  filter.  Matched  to  a  predetermined 
section  of  the  received  signal,  the  matched  filter  outputs 
periodic  amplitude  peaks  indicating  the  occurrence  of  a 
subsequence  of  the  desired  code.  Kt  this  moment  the  local  PN 
sequence  generator  is  started  from  the  corresponding  initial 
state.  At  the  same  instant  the  tracking  loop  is  activated 
(12:5.9.2) . 

Matched  filters  in  this  application  typically  employ 
surface  acoustic  wave  (SAW)  technology.  SAW  devices  are 
ideal  choices  for  use  in  spread  spectrum  communication 
systems  which  must  process  high  data  rate  signals.  The 
SAW*>tapped  delay  line  matched  filter  is  a  passive  correlation 
technique  which  maximizes  the  filter  output  peak  signal-to- 
mean  noise  ratio.  However,  the  length  of  SAW  matched  filters 
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is  United  by  propagation  losses  and  the  filter  cannot  be 
Inplenented  to  match  the  length  of  typical  PN  code  sequences 
(13:230) . 

The  system  of  Interest  presents  a  combined  matched 
fllter/serlal  search  rapid  acquisition  strategy  In  a  direct 
sequence,  multiple-access  environment.  A  commonly  used 
method  for  PN  code  acquisition  In  direct  sequence  spread 
spectxrum  systems  Is  the  serial  search  of  the  relative 
positions  of  the  received  signal  and  the  locally  generated 
reference  waveform  in  discrete  steps  of  a  fraction  of  a  code 
chip.  If  no  a  priori  timing  Information  Is  available  to  the 
receiver,  uniform  search  over  the  long  pseudonoise  code 
periods  lead  to  acquisition  times  that  are  prohibitive. 

Prior  information  about  the  approximate  position  of  the 
received  code  leads  to  a  decrease  in  the  mean  acquisition 
time  by  several  orders  of  magnitude  over  a  vide  range  of  the 
receiver  input  signal-to-noise  ratios  (14:157).  In  this 
system,  the  initial  uncertainty  is  reduced  with  matched 
filter  detection  of  shorter  code  segments  of  the  received  PN 
sequence.  Attention  is  focused  on  initial  acquisition.  In 
other  words,  the  initial  establishment  of  timing 
synchronization  from  total  timing  uncertainty. 

Signal  Elan  And  Qp?rati9n 

The  most  straightforward  acquisition  strategy  for  long 
pseudonoise  sequences  uses  active  correlation  and  a  simple 
serial  search  algorithm  to  determine  the  correct  code  epoch. 
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This  Bsthod  is  ths  standard  for  comparison  of  a  system  that 
augments  this  technique  with  an  auxiliary  matched  filter 
channel  to  speed  code  acquisition.  A  block  diagram  of  the 
receiver  synchronization  section  is  shown  in  Figure  4.  The 
major  components  of  the  system  are  the  pseudonoise  sequence 
generator,  the  matched  filter  channel,  and  a  quadrature 
correlator  as  the  active  correlation  device. 


Figure  4.  Synchronization  System 


The  matched  filter  channel  shown  in  Figure  5  consists  of 
the  matched  filter,  envelope  detector,  accumulator  device, 
and  a  threshold  detector.  The  matched  filter  can  be 
implemented  with  the  aid  of  tapped  delay  lines.  The  most 
suitable  IF  devices  for  use  in  tapped  delay  line  correlators 
tend  to  be  the  surface  acoustic  wave  type  as  discussed 
earlier.  Figure  6  shows  a  schematic  of  a  tapped  delay  line 
for  use  as  a  matched  filter  in  spread  spectrum  receivers-a 
usable  filter  would  contain  many  delay  elements. 
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The  sequence  PN(t)  is  a  saxinal  length  sequence  whose 

autocorrelation  function  R^(t)  at  |r|  2  is  very  low. 

The  chip  period  T_  is  selected  in  such  a  way  that  it  is  an 

c 

integer  sultiple  of  the  carrier  period.  The  delay  line  has  M 
taps  spaced  at  tine  intervals  of  T^.  The  delay  intervals 
must  accurately  represent  the  clock  period  of  the  code 
sequence  to  be  detected  (12:5.9.2). 

The  taps  of  the  delay  line  matched  filter  are  programmed 
such  that  the  matched  filter  will  recognize  a  particular 
svibsequence  of  length  M  of  the  incoming  PN  sequence.  The  tap 
values  will  pass  the  signal  unaltered  or  with  a  phase 
reversal  in  accordance  with  the  desired  section  of  sequence 
PN(t).  With  a  length  of  M  chips,  k^^  of  the  M  output  signals 
of  the  delay  line  are  fed  directly  to  the  adding  network  and 
the  remainder  k^  «  M  •  )c^  after  a  phase  reversal  of  180*. 
Most  of  the  time,  assuming  the  PN  code  sequence  signal  levels 
are  equally  probable,  the  amplitude  of  the  summation  is  on 
the  order  of  JlP.  When  the  phase  reversals  of  the  input 
signal  match  those  reversals  programmed  into  the  delay  line, 
the  output  summation  is  on  the  order  of  N VJP  since  M  signals 
having  the  same  phase  are  added  (12:5.9.2). 

The  rectified  and  filtered  periodic  peaks  at  the  output 
of  the  envelope  detector  mark  a  certain  state  of  the  received 
PN  sequence.  This  signal  is  therefore  suitable  as  a 
synchronization  pulse  for  setting  and  starting  of  the  local 
PN  generator.  There  is,  however,  more  to  consider. 
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The  subsequence  of  PN(t)  progranmed  into  the  tapped 
delay  line  nust  be  chosen  such  that  the  sain  correlation  peak 
exceeds  the  partial  autocorrelation  characteristics.  In 
practical  systems,  the  input  signal  to  the  matched  filter 
contains  a  narrow-band  data  modulation.  This  impairs  the 
operation  of  the  matched  filter.  The  impairment  can  be  made 
negligibly  small  as  long  as  the  bandwidth  of  the  data 
modulation  is  relatively  narrow.  The  signal-to-noise  ratio 
(SNR)  of  an  interfering  signal  at  the  output  of  the  delay 
line  is  higher  by  about  a  factor  M  as  compared  to  its  input 
SNR.  The  setting  and  starting  accuracy  is  impaired  by  this 
interference.  The  Interference  can  offset  synchronization 
epochs  and  prevent  detection  of  correct  synchronization 
pulses.  The  rate  of  false  alarms  and  loss  of  correct 
synchronization  pulses  become  greater  as  the  input  SNR  of  the 
target  transmitter  becomes  smaller  and  as  the  input  deviates 
from  the  desired  code  due  to  data  modulation  (12:5.9.3). 

In  order  to  ensure  acceptable  values  in  the  probability 
of  detection  (P^)  and  the  probability  of  false  alarm  (P^jj) » 
the  accumulator  circuit  processes  a  designated  number  of  main 
correlation  peaks  prior  to  threshold  detection  and  initiation 
of  the  PN  generator.  The  number  of  correlation  peaks 
processed  is  based  on  the  peak  SNR  at  the  output  of  the 
matched  filter.  The  threshold  detect  circuit  is  set  to 
initiate  the  PN  generator  such  that  the  active  correlator 
will  synchronize  in  the  first  dwell  period  for  given 
probabilities  of  false  alarm  and  detection. 
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Analytical  Evaluation 


The  analysis  of  the  quadrature  correlator  and  the 
matched  filter  channel  serves  two  purposes:  1)  the  output 
signals  must  be  characterized  in  order  to  verify  the 
simulation  models,  and  2)  the  mean  acquisition  time 
performance  must  be  assessed  to  gauge  the  accuracy  of  the 
simulation  results. 

The  first  objective  is  to  establish  the  performance  of 
the  quadrature  correlator.  To  assess  the  improvement  in 
synchronization  time  provided  by  the  auxiliary  matched  filter 
channel,  the  performance  of  the  quadrature  correlator  as  a 
stand-alone  synchronization  circuit  is  analyzed.  This  will 
serve  as  the  standard  to  judge  improvement  in  acquisition 
time  via  the  auxiliary  matched  filter  channel.  The 
performance  due  to  inclusion  of  the  matched  filter  channel 
will  then  be  determined. 

Quadrature  Correlator.  The  quadrature  correlator  shown 
in  Figure  7  is  a  coherent  single  dwell  acquisition  system. 

It  is  assumed  that  the  receiver  is  capable  of  determining 
good  estimates  of  the  carrier  phase  produced  by  the 
propagation  delay  and  Doppler.  The  received  signal  plus 
noise  is  demodulated  by  in-phase  and  quadrature  carrier 
reference  signals  whose  frequency  is  equal  to  that  of  the 
transmitter  carrier  frequency  but  with  an  arbitrary  phase. 

The  demodulated  and  despread  signals  are  passed  through 
integrators  with  identical  dwell  times.  It  is  assumed  that 
the  rectifiers  are  square-law  devices  to  remove  the  unknown 


24 


Figure  7.  Quadrature  Correlator 


data  modulation  and  unknown  carrier  phase  (11:67).  It  is 
also  assumed  that  the  output  of  the  correlator  is  fed  to  a 
simple  threshold  device  controlling  an  implementation  of  a 
serial  search  strategy.  The  evaluation  will  proceed  with  a 
characterization  of  the  quadrature  correlator  output  signal. 

Correlator  Output  Signal.  The  received  modulated 
signal,  in  the  absence  of  noise,  is  of  the  form 

r(t)  -  y2Pdj(t)PNj(t)cosw^t  +  y2PdQ(t)PNQ(t)sinw^t  (12) 

where  d^(t)  and  dQ(t)  are  the  in-phase  (I)  and  quadrature  (Q) 
data  streams,  and  PN^(t)  and  PNQ(t)  are  the  pseudonoise 
sequences  with  the  I  signals  on  a  carrier  shifted  90*  from 
the  Q  carrier  (15:347). 

A  -conventional  QPSK  signal  is  assummed.  The  PN  code 


epoch  times  are  common  to  both  the  I  and  Q  signals  and 


without  indication  otherwise  in  the  Snirnov  description,  the 
I  and  Q  pseudonoise  sequences  are  assumed  identical.  The 
receiver's  replica  signals  are: 

Sj(t)  -  y2PPNj(t+T)cos(w^t+0)  (13) 

SQ(t)  -  y2PPNQ(t+r)sin(w^t+«)  (14) 

where  r  is  the  time  offset  of  the  PN  replica  sequence  and  0 
is  the  carrier  reference  offset.  The  received  signals  are 
modulated  by  the  replica  signals  and  input  to  the  integrate 
and  dump  circuits.  The  input  to  the  in-phase  integrator  is 


Xj(t)  *  r(t) •Sj(t) 

«  2P  j^dj  (t)  PNj  (t)  PNj  (t+T )  cosw^t  •  cos  (v^t+<f>) 

+  dQ(t)PNQ(t)PNj(t+r)sinw^tcos(w^t+0)j  (15) 

Given  the  period  of  integration  is  over  a  complete  PN 
sequence  and  assiuning  the  data  rate  is  low  enough  relative  to 
the  chip  rate  that  dj(t)  and  dQ(t)  are  constants  over  this 
period,  the  output  of  the  in-phase  integrator  channel  is 


rnTc 

Xj(t)dt 


rnT. 


p  dj(t) 


PNj (t) PNj(t+T)COS0[l+COS2W^t)dt 


+  P  <jQ(t) 


.nT^ 

PNq (t) PNj (t+f ) sin0[ l-cos2WQt )dt 
0 


(16) 


Since  the  period  of  integration  is  an  integral  multiple  of 
the  period  of  the  carrier,  the  integrate  and  dump  circuit 
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acts  as  a  low  pass  filter.  The  result  is  as  follows: 


f 


nTc 

Xj(t)dt 

0 


Slmlliarly, 


rnTj, 

XQ(t)dt 


-  P  nT^  Rpjj(r)  j^dj(t) -cos^  +  dQ(t)  sin0j 
for  the  quadrature  channel: 

-  P  nT^  Rpjj(T)  [^dj(t)  sin^  +  dQ(t)  cos0j 


(17) 


(18) 


Assuming  the  rectifiers  are  square-law  devices,  the  data 
modulation  and  the  unknown  carrier  phase  is  effectively 
removed  and  the  output  of  the  correlator  is: 


Z^(t)  .  2[p  nT^  Rp„(f)] 
With  the  bit  energy  as  unity, 
Z^(t)  .  2[nRp„(r)] 


(19) 


(20) 


When  r  equals  zero,  the  reference  PN  sequence  is  synchronized 

with  the  received  PN  sequence  and  Z  (t)  achieves  its  maximum 

c 

value;  in  other  words,  a  main  correlation  peak  occurs  at  the 
output  of  the  correlator. 

g9xrslat9c  Acqui8iti9n  Zise  P?rff9]rwan99.  it  is 
assumed  that  the  correlator  output  signal  is  applied  to  a 
simple  threshold  device  which  in  turn  drives  a  serial  search 
strategy.  If  the  correlator  output  is  above  the  threshold 
level,  then  the  search  algorithm  enters  the  code  tracking 
loop  mode.  The  threshold  indication  is  verified  in  the 
tracking  mode  with  an  extended  dwell  time.  The  extended 
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dwell  time  is  an  integral  multiple  of  the  dwell  time,  where 
the  dwell  time  is  the  correlator  period  of  integration.  If 
the  verification  is  successful,  then  the  true  code  phase 
position  has  been  acquired  and  the  search  comes  to  an  end. 

If  not,  the  search  continues  and  the  time  penalty  for  the 
false  alarm  is  the  extended  dwell  time  of  the  verification 
process.  If  the  correlator  output  is  below  the  threshold 
value,  the  local  PN  generator  is  stepped  to  its  next  position 
and  the  search  proceeds  (11:17). 

It  is  typical  practice  to  require  the  alignment  between 
the  received  and  local  PN  codes  be  within  one-half  of  a  code 
chip  period  before  entering  the  tracking  mode.  The  PN  code 
generator  is  stepped  in  one-half  chip  periods  over  the 
uncertainty  region.  In  the  absence  of  a  priori  information 
concerning  the  code  phase  position,  the  uncertainty  region 
between  the  received  PN  code  and  the  local  PN  code  could  be 
as  much  as  the  entire  code  period  (11:18). 

Smirnov  defines  the  region  of  uncertainty  as  the  entire 
code  period.  Smirnov  presents  Eq  (21)  as  the  time  for  entry 
into  synchronization  for  the  quadrature  correlator  (3:68). 

V  “  '21) 

where  n  is  the  number  of  chips  in  a  full  PN  code  sequence, 

T  is  the  chip  duration,  and  E  is  the  number  of  frequency 
intervals  in  the  region  of  uncertainty.  This  result  is  taken 
as  the  average  time  for  acquisition  for  complete  timing 
uncertainty  and  the  simple  case  where  probability  of 
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detection  is  unity,  probability  of  false  alam  is  zero,  and 
the  dwell  time  (as  it  turns  out)  is  for  full  period 
correlation . 

This  result  can  be  verified  using  the  aean  acquisition 
time  of  a  serial  direct  sequence  search  system  (11:22).  The 
expression  given  below  assumes  the  search  is  performed  in 
half-chip  increments  (cells)  and  in  the  absence  of  code 
Doppler  and  carrier  Instabilities. 


acq 


(2  -  P^)(l  K  P^^) 
2P^ 


q.r 


(22) 


As  stated  above,  the  case  given  by  Smirnov  sets  the 
probability  of  detection  (P^)  to  unity  and  the  probability  of 
false  alarm  (P^^)  to  zero.  With  the  dwell  time  set  for  full 
period  correlation,  the  number  of  cells  to  search  is  q  «  2n 
where  n  «  M  ,  the  total  number  of  chips  in  a  full  period. 

It  follows  that  the  integration  time  for  a  full  period  is 
the  result  of  Eq  (21)  is  repeated: 


T 


acq 


n^T. 


(23) 


This  result  is  the  mean  acquisition  time  for  searching  the 
frequency  interval  in  which  the  code  is  located  (E  1)  . 

Before  proceeding  with  the  analysis  of  the  matched 
filter  channel,  the  standards  for  correlator  performance  must 
be  expressed  in  terms  of  receiver  operating  characteristics. 
When  the  matched  filter  Initiates  the  PN  generator,  it  does 
so  such  that  for  given  probabilities  of  false  alarm  and 
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detection,  the  quadrature  correlator  will  synchronize  upon 
detection  of  a  single  pulse  in  the  first  dwell  period.  The 
ability  of  the  quadrature  correlator  to  successfully  indicate 
synchronization,  as  specified  by  the  above  probabilities,  is 
a  function  of  its  output  SNR.  This  data  is  tabulated  for 


single  and  multiple  pulse  detection  of  signals  and  will  be  ^ 

the  standard  to  determine  the  required  matched  filter  output 
SNR  (16:205, 248>261) .  If  the  matched  filter  initiates  the  PN 
generator  with  the  specified  output  SNR,  the  correlator  will  I 

synchronize  with  the  same  probabilities. 

Figure  8  specifies  average  output  SNR  requirements  for 
single  pulse  detection  over  a  wide  range  of  and  P^.  All  .  .Hl 

further  analysis  will  develop  with  the  peak  SNR  rather  than 
the  average  SNR.  The  values  provided  will  be  scaled  by  -t-SdB. 
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Matched  Filter  Channel .  The  main  advantage  of  natched 
filtering  in  a  synchronization  circuit  is  the  potential  of 
very  fast  acquisition  times.  The  IF  matched  filter  placed  in 
an  auxiliary  channel  continuously  correlates  a  subsequence  of 
the  incoming  signal  with  that  of  the  reference  subsequence  of 

► 

the  filter.  A  correlation  peak  appears  at  the  output  of  the 
matched  filter  when  the  epoch  of  the  incoming  signal  matches 

with  that  of  the  filter.  The  timing  signal  is  obtained  by 

I 

rectifying  the  matched  filter  output  and  accumulating  a 
number  of  successive  correlation  peaks  until  a  threshold 
level  is  crossed.  The  timing  signal,  initiated  at  the 

» 

instant  of  the  threshold  crossing,  sets  the  epoch  of  the 
local  PN  generator  and  receiver  synchronization  is  achieved. 

Matched  Filter  Q.vi.tPUt  Signal.  As  given  in  Eq  (12)  , 

I 

the  received  modulated  signal,  in  the  absence  of  noise,  is  of 
the  form: 

r(t)  *  y^  dj(t)PNj(t)cosw^t  +  y^  dgCt) PN^Ct) sinw^t  (12) 

where  d^(t)  and  dQ(t)  are  the  I  and  Q  data  streams,  and 

PN^(t)  and  PNg(t)  are  identical  pseudonoise  sequences  with 

common  epoch  times  for  the  I  and  Q  channels.  The  pseudonoise 

signals  have  a  clock  frequency  f  *  1/T  with  a  period 

c  c 

equal  to  n-T^  where  n  is  the  total  number  of  chips  in  a  full 
sequence . 

Assuming  additive  white  Gaussian  noise  at  the  input  to 
the  receiver,  the  impulse  response  of  the  matched  filter  is 
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h(t) 


y2P-  PN-.  (MT^-t)  cosw  t 
X  c  o 

0 


0  ^  t  ^  MT^ 
c 

otherwise 


(24) 


The  output  of  the  matched  filter  is  the  convolution  of  the 
received  signal  and  the  impulse  response  of  the  filter  as 
given  by  u(t)  »  r(t)  *  h(t)  .  Assuming  the  data  rate  is  low 
relative  to  the  chip  rate,  dj(t)  and  dgCt)  are  constants  over 
the  period  of  integration  for  the  convolution  and  the  output 
is: 


u(t) 


rMT 


c 


Pdj(t) 


PNj ( T ) PMj [ r+ ( t-MT^) ] cosw^t [ 1+cos2w^t ] dr 

Jq 


rMT 


c 


+  p  dQ(t) 


PNq ( r ) PN j [ T+ ( t-MT^) ] sinw^t [ l-sin2w^T ] dr 
0 


(25) 


Assuming  the  period  of  integration  is  much  greater  than  the 
period  of  the  carrier,  u(t)  is  given  by 


u(t)  »  P  MT^  Rpjj(t-MT^)  [dj(t)cosw^t  +  dQ(t)sinw^t]  (26) 


The  output  of  the  matched  filter  is  passed  through  an 
envelope  detector  consisting  of  a  square-lav  device  and  a  low 
pass  filter.  The  data  modulation  is  effectively  removed 
through  the  square-law  device  with  the  following  output: 


Cu(t)]2 


(27) 
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With  the  bit  energy  as  unity,  the  input  to  the  threshold 
detector  is  given  by 


LPF[u*(t) ] 


|t|  <  MT^ 
|t|  >  MT^, 


(28) 


When  |t|  ^  MT^  ,  the  reference  PN  sequence  is  synchronized 
with  the  received  PN  sequence  and  the  output  of  the  matched 
filter  channel  achieves  its  maximum  value,  in  other  words,  a 
full  correlation  peak  occurs.  At  all  other  times,  n^(t) 
represents  the  output  of  partial  autocorrelation  peaks  that 
occur  due  to  the  correlation  characteristics  of  subsequences 
of  the  selected  PN  code.  This  value  cannot  be  determined 
analytically  but  will  be  quantified  via  simulation. 

Watghgd  Ultsi:  Acquisition  Tims  pgrfprB!anse.t  m 
addition  to  the  main  correlation  peak  and  the  partial 
autocorrelation  peaks  due  to  the  PN  synchronization  waveform, 
additional  signals  appear  at  the  output  of  the  detector.  The 
IF  bandwidth  includes  thermal  noise,  the  data  waveform  of  the 
desired  signal,  and  in  a  multiple-access  environment-the  data 
and  PN  code  waveforms  of  other  users.  These  additional 
signals  contribute  to  the  partial  autocorrelation  output  and 
collectively  degrade  the  output  signal-to-noise  ratio. 

As  stated  previously,  the  driving  performance 
characteristic  is  the  matched  filter  peak  output  signal-to- 
noise  ratio  •  1'°  determine  the  average  signal 

power  for  each  of  the  above  interfering  signals  and  the  peak 
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I 


power  of  the  desired  PN  code  sequence  must  be  determined. 

The  performance  of  the  matched  filter  channel  can  be 
characterized  in  terms  of  the  first  two  moments  of  the 
probability  density  function  for  the  amplitude  of  the  matched 
filter  output.  The  partial  correlation  output  of  the  matched 
filter  is  treated  as  a  Gaussian  random  variable  at  the  input 
to  the  square-law  detector.  The  basis  for  the  Gaussian 
assumption  is  the  relationship  between  the  matched  filter 
length  and  the  overall  PN  code  length.  Because  the 
correlation  time  of  the  matched  filter  (MT_)  is  much  less 
than  the  PN  code  time  period  (n-T  ),  the  correlation  output 

w 

amplitude  has  a  binomial  distribution.  With  the  additional 
relationship  that  the  correlation  period  is  much  greater  than 
unity,  the  binomial  distribution  behaves  like  a  Gaussian 
distribution  (11:78). 

The  PN  code  correlation  function  is  triangular  over  an 
interval  of  plus  and  minus  one  code  chip  (-T  ,T  ) .  Let  H. 

w  C  X 

denote  the  hypothesis  that  the  received  code  waveform  and  the 
reference  waveform  of  the  matched  filter  are  misaligned  by 
less  than  one  code  chip.  The  alternate  hypothesis 
denotes  that  the  relative  alignment  is  greater  than  or  equal 
to  a  code  chip  period  (11:61). 

The  input  chip  waveform,  y^  r(t+T),  correlated  by  the 
matched  filter  reference,  c(t+T')  over  the  period  MT  is 
given  by  (11:77): 


C 


pMTc 

r(t+T)  •c(t+T')dt 


(29) 
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The  expected  value  of  the  correlation  conditioned  upon 
whether  the  incoming  chip  vavefom  is  correlated  with  the 
reference  waveform  (H^)  or  uncorrelated  (H^)  is  given  by 
(11:77) : 


e{c1h^} 


ySp-HT^- (i-|n|) 

0 


i  -  1 
i  -  0 


(30) 


where  n  is  the  fractional  timing  offset  with  respect  to  a 
chip  between  the  two  code  waveforms  as  given  by, 

±  (31) 

and  is  the  smallest  integer  such  that  n  lies  in  the 
interval  (-1,-fl).  It  is  apparent  that  the  expected  value 
reaches  a  peak  value  of  y^  MT_  when  the  codes  are  correlated 
(n=0) .  since  partial  correlation  is  a  Gaussian  random 
variable  with  zero  mean,  the  expected  value  of  uncorrelated 
code  segments  is  zero. 

The  power  of  the  uncorrelated  waveforms  is  expressed  by 
the  variance  given  below  (11:77): 

var(C|H. )  «  2PMT^^G. (n)  (32) 

X  w  X 


where 


G.  (0) 


n 


2 


1  -  2|n|  +  20^ 


i  »  1 
i  «  0 


(33) 
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other  users  will  occur  with  an  offset  outside  the  correlation 


interval.  As  seen  in  Figure  9,  the  variance  for  an 
uncorrelated  vavefom  repeats  every  chip  period  outside  the 
correlation  interval.  Assuming  that  the  timing  offset  for 
these  waveforms  is  uniformly  distributed  over  a  chip  period, 
the  average  offset  will  be  n  «  2/3  as  shown  below: 


n 

I 


pi 

(l-2|a|+2a^)do  «  2/3 
**0 


(34) 


The  average  power  for  the  data  waveforms  associated  with 
other  users  is  therefore  2/3-P^,MT^  .  The  average  power 
for  the  PN  code  waveforms  associated  with  other  users  is 

2/3P^,MT^2^ 

With  thermal  noise  at  the  input  to  the  matched  filter 

expressed  as  N  «  N^BW^p  and  no  correlation  with  the 

reference  signal  (n»l  or  -1) ,  the  noise  power  at  the 

2 

output  of  the  matched  filter  channel  is  N-MT  .  The  total 
average  power  of  Interfering  sources  can  now  be  expressed: 


tot 


"^c^  2/3  K(Pjj,  +  P^,)] 


(35) 


Where  K  ^  0  is  the  number  of  interfering  users. 

The  peak  signal  power  of  the  desired  PN  code  waveform  is 
the  square  of  the  peak  (n«0)  expected  value  for  a  correlated 
waveform  and  is  equal  to  2P^{MT^)^.  The  peak  SNR  at  the 
output  of  the  matched  filter  channel  can  now  be  expressed: 
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(36) 


2P-M 

X - C - 

^  N  +  Pjj  +  2/3K(Pjj,  +  P^,) 

This  signal-to-nolse  ratio  will  detennina  the  nimber  of 
correlation  peaks  that  must  be  processed  to  achieve 
synchronization  with  given  probabilities  of  false  alarm  and 
detection.  To  achieve  synchronization  in  the  first  code 
period,  the  matched  filter  ouput  signal-to-noise  ratio  must 
meet  the  same  single-pulse  detection  requirements  given  in 
Figure  8  for  the  quadrature  correlator.  For  a  probability  of 
detection  equal  to  0.9  and  selected  values  for  probability  of 
false  alarm,  the  peak  output  signal-to-noise  ratio 
requirements  are  summarized  below  in  Table  1. 


Table  1.  Single-Pulse  Detection  Requirements 


Pfa 

Peak  Output  SNR 
Requirements 

lo"^ 

13.7  dB 

10”5 

15.5  dB 

o 

1 

16.2  dB 

As  evident  in  Eq  (36),  the  matched  filter  output  signal- 
to-noise  ratio  is  dependent  on  the  length  of  the  filter,  the 
number  of  users  in  a  multiple-access  system,  and  the  level  of 
thermal  noise.  Given  these  factors,  the  matched  filter  may 
not  achieve  the  required  output  signal-to-noise  ratios 
required  for  single-pulse  detection.  At  lower  output  signal- 
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to-nolse  ratios,  several  correlation  spikes  Bust  be  processed 
prior  to  threshold  detection  and  initiation  of  the  local  PN 
generator.  Requireaents  for  multiple  pulse  detection  of 
signals  are  tabulated  (16:248-261)  and  are  summarized  below. 
Figure  10  gives  the  number  of  pulses  required  for  detection 
at  specified  probabilities  of  false  alarm  and  detection  over 
a  range  of  peak  output  slgnal-to-noise  ratios. 


Ntroer  of  Pulses  for  Po  =  o  3 
O  l0e-3  +  l0e-6 


Figure  10.  Multiple-Pulse  Detection 
Requirements  (16:246-261) 

The  threshold  level  at  the  output  of  the  accumulator 
controls  the  number  of  pulses  processed  by  the  matched  filter 
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prior  to  the  initiation  of  the  local  PN  generator.  At  the 
end  of  each  PN  code  period,  the  natched  filter  outputs  a  main 
correlation  peak  of  amplitude  on  the  order  of  the  length  of 
the  filter.  This  main  correlation  spike  and  the  associated 
noise  throughout  the  remaining  PN  code  period  is  squared 
through  the  envelope  detector  and  input  to  the  accumulator. 
Dependent  on  the  noise  level  accompanying  the  main 
correlation  peak,  the  threshold  level  is  set  to  require  the 
specified  number  of  pulses  to  accumulate  prior  to  threshold 
crossing.  If  the  threshold  requirements  are  met,  the  local 
PN  generator  will  be  initiated  such  that  the  quadrature 
correlator  will  indicate  synchronization  in  the  first  dwell 
period  with  corresponding  probabilties  of  false  alarm  and 
detection. 

Acquisition  Zims  Performance 

The  matched  filter  output  SNR  and  therefore  the 
acquisition  time  of  the  matched  filter  channel  is  influenced 
by  the  length  of  the  filter,  the  number  of  interfering  users, 
and  the  level  of  thermal  noise,  in  this  section,  each  of 
these  influences  will  be  examined  in  order  to  characterize 
the  acquisition  time  performance  of  the  system.  A  comparison 
will  then  be  performed  to  highlight  the  improvement  in 
acquisition  time  offered  by  the  auxiliary  matched  filter 
channel  against  the  quadrature  correlator  alone. 

The  first  step  is  to  simplify  the  equation  describing 
the  peak  output  SNR  of  the  matched  filter.  In  Eq  (36) ,  the 
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power  of  the  data  and  PN  code  wavefoms  will  all  be  set  equal 


such  that  -  Pjj  ■  Pg,  *  Pjjt  •  Allowing  the  combined  power 
of  the  data  and  PN  code  waveforms  of  a  particular  transmitted 
signal  to  be  represented  as  S  «  P^  -t-  P^  ,  the  matched  filter 
peak  output  SNR  reduces  to: 

M«i 

i -  (37) 

1  +  S^(l/2  +  2/3  K) 

Where  *  S/N  is  the  matched  filter  input  SNR  of  a  single 
transmitted  signal  and  K  is  the  nvunber  of  other  users. 

As  stated  previously,  for  output  signal-to-noise  ratios 
less  than  the  requirements  for  single-pulse  detection,  the 
matched  filter  channel  must  process  a  number  of  pulses  to 
meet  a  particular  false  alarm  rate  and  probability  of 
detection.  Recognizing  that  each  pulse  the  matched  filter 
channel  processes  equates  to  a  PN  code  period,  the 
requirements  for  multiple-pulse  detection  given  in  Figure  10 
can  be  used  to  evaluate  acquisition  time  performance. 

Acquisition  time  performance  in  terms  of  code  periods  is 
specified  by  the  number  of  pulses  processed  prior  to 
initiation  of  the  local  PN  generator  and  the  length  of  code. 
The  length  of  a  code  period  is  equal  to  the  code  chip  time 
multiplied  by  the  number  of  chips  in  the  code  sequence  or 
n-T^.  Letting  N^  represent  the  number  of  pulses  processed 
according  to  the  multiple-pulse  detection  requirements  of 
Figure  10,  the  mean  acquisition  time  is  given  by: 


'MA 


p  C 


(38) 
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Acquisition  tine  performance  for  a  range  of  input 
signal-to-noise  ratios,  matched  filter  lengths,  and  number  of 
other  users  is  given  below.  Figures  11  through  15  show 
performance  characteristics  for  an  ideal  false  alarm  rate 
■  lo"®  and  detection  probability  P(j  *  0*®  while 

Figures  16  through  20  show  performance  characteristics  at 

“3  "3 

»  10  .  A  probability  of  false  alarm  at  10  can  be 

simulated  in  reasonable  time  frames  and  the  expected  values 

in  Figures  16  through  20  are  provided  to  highlight  the 

performance  characteristics  that  will  be  simulated. 

Improvement  in  acquisition  time  offered  by  the  matched 

filter  channel  is  evident  in  the  ratio  of  the  mean 

acquisition  times  for  the  correlator  and  that  of  the  matched 

filter  given  in  Eg  (39).  Recalling  Eg  (22)  and  recognizing 

that  KPj^  «  1  for  the  case  of  interest  (P^^  •  lo“^) ,  the 

mean  acquisition  time  for  the  correlator  acting  alone  is: 

-  1-22  (39) 

where  n  is  the  number  of  code  chips  in  the  PN  code  period. 

The  acquisition  time  is  improved  by  a  factor  (I)  expressed  in 
the  following  ratio: 

n 

I  -  1.22 -  (40) 

"p 

As  shown  in  Figures  11  through  15,  for  signal-to-noise 
ratios  in  the  range  above  0  dB,  a  variation  of  the  input  SNR 
has  no  effect  on  the  acquisition  times.  In  this  range,  the 
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input  SNR 


Ir^ut  SNfi  CdBi  1  »1)  I  input  COB) 
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Figure  19.  MF  Length  =  512  »  lO’^/P^  =  0.9) 


Acquisition  Tm»  (»  Code  Perioosj 
O  t  *  0  .KiA  OC.tE  kJB«rs 


fa  -  -  o-®* 


Figure  20.  MF  Length  »  1024  (P 


primary  Influence  on  acquisition  time  performance  for  a  given 
matched  filter  length  is  the  niunber  of  users  in  the  system. 
Improvement  in  the  acquisition  times  over  this  range  of  input 
signal-to-noise  ratios  is  sho%ni  in  Figure  21  for  a  code 
length  of  n  *  2047  chips. 


As  the  input  signal-to-noise  ratio  decreases  below  0  dB, 
the  number  of  users  in  the  system  becomes  less  of  an 
influence  on  the  acquisition  times.  For  further  decreases  in 
the  input  signal-to-noise  ratio,  the  curves  presented  in 
Figures  11  through  15  asymptotically  approach  a  doubling  in 
the  acquisition  time  for  each  1.5  dB  decrease  in  the  input 
SNR  (16:262).  For  the  input  SNR  range  below  -20  dB,  noise  is 
the  dominant  influence  on  acquisition  time.  To  demonstrate 
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this  trend,  the  Improvement  In  acquisition  time  for  an  input 
SNR  of  -10  dB  is  given  in  Figure  22. 
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Figure  22.  Acquisition  Time  Improvement  (Input  SNR  »  -10  dB) 

As  evident  in  Figure  22,  increasing  the  number  of  users 
does  not  produce  the  degree  of  variation  in  the  improvement 
factor  that  is  seen  at  higher  input  signal-to-noise  ratios. 

In  addition,  the  improvement  levels  are  significantly  reduced 
to  due  the  high  AHGN  noise  levels. 

To  summarize,  it  is  the  peak  output  SNR  of  the  matched 
filter  that  determines  acquisition  time  performance.  A 
number  of  factors  influence  this  system  parameter.  First  of 
all,  the  partial  autocorrelation  of  the  PN  code  of  interest 
outputs  correlation  noise.  Other  influences  on  the  output 
SNR  include  the  cross-correlations  of  the  data  waveform  of 


49 


the  signal  of  Interest,  the  data  and  PM  code  waveforms  of  the 
other  users,  and  thermal  noise.  The  amplitudes  of  these 
signals  are  assumed  to  have  a  Gaussian  distribution.  The 
influences  of  the  Gaussian  signals  and  the  length  of  the 
matched  filter  are  quantified  in  Eg  (36)  for  the  matched 
filter  peak  output  SNR. 

For  a  selected  and  P^,  the  peak  output  SNR 
determines  whether  a  signal  can  be  detected  within  the 
tabulated  single-pulse  or  multiple-pulse  detection 
requirements.  If  the  output  SNR  does  not  meet  single-pulse 
detection  requiremets,  the  number  of  pulses  that  must  be 
processed  to  meet  the  multiple  pulse  detection  requirements 
determines  the  acquisition  time.  The  time  between  successive 
pulses  is  equal  to  a  full  PN  code  period.  Therefore,  the 
acquisition  time  performance  is  directly  related  to  the 
matched  filter  output  SNR  as  a  function  of  the  matched  filter 
length,  the  number  of  other  users  in  the  system,  and  the  AWGN 
level . 

The  next  chapter  proceeds  with  modeling  the  system 
components  using  the  Block-Oriented  Systems  Simulator.  The 
output  signal  analysis  in  this  chapter  will  be  used  to 
validate  the  simulation  models  for  the  matched  filter  channel 
and  the  quadrature  correlator.  The  single  and  multiple-pulse 
detection  requirements  will  be  used  to  select  simulation 
times  and  associated  threshold  detection  levels  to  test 
performance  of  the  matched  filter  via  simulation. 
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I2L.  BOSS  Mpleling  stt  Syst^w 


Th«  Block  Oriented  Systems  Simulator  (BOSS)  provides  for 
simulation-based  analysis  and  design  of  any  system  which  can 
be  represented  in  block  diagram  form.  Simulation  models  of 
systems  are  constructed  using  a  heirarchial  block  diagram 
approach  where  each  block  performs  a  signal  processing 
operation.  Once  a  correct  description  of  the  system  is 
expressed  in  block  diagram  form,  BOSS  supports  a  time-domain 
(waveform  level)  simulation  (1:36.1.1).  The  pseudonoise  (PN) 
synchronization  system  to  be  simulated  is  composed  of  a 
single  dwell  quadrature  correlator  augmented  by  a  matched 
filter  channel  to  speed  PN  code  acquisition.  The  simulation 
model  will  be  configured  to  analyze  the  initial  acquisition 
characteristics  of  the  matched  filter  channel. 

The  chapter  begins  with  the  approach  to  modeling  the 
synchronization  system.  The  second  section  is  dedicated  to 
the  successive  steps  and  analysis  to  perform  a  bottom-up 
modeling  of  the  individual  subsystem  components  using  BOSS 
functional  modules.  The  final  section  will  discuss 
integration  of  the  simulation  models  for  a  single  user  system 
and  a  multiple-access  system. 

Top-down  Decomposition 

The  first  step  in  the  BOSS  modeling  process  is  to 
perform  a  top-do%m  decomposition  of  the  system.  The  major 
components  of  the  system,  the  (1)  Pseudonoise  Sequence 
Generator,  the  (2)  Matched  Filter  Channel,  and  the 


(3)  Quadrature  Correlator  consist  of  the  block  diagram 
representations  shown  below  in  Figure  23. 


Figure  23.  Synchronization  System  (3:68) 

Further  decomposition  of  the  functional  blocks  is  necessary 
to  suppport  BOSS'S  bottom-up  hierarchial  concept  of  system 
design.  BOSS  supplies  a  library  of  low-level  modules  called 
primitives  to  construct  more  complex  modules.  Nodules  of 
increasing  complexity  can  in  turn  be  used  in  higher  level 
modules  in  a  bottom-up  hierarchial  fashion  until  the  visual 
representation  of  the  simulation  matches  the  high-level 
system  block  diagram  of  Figure  23. 

The  first  simulation  model  will  be  the  PN  sequence 
generator.  The  generator  will  then  be  used  in  the 
construction  of  the  QPSK  direct  sequence  transmitter.  Once 
the  transmitter  has  been  successfully  tested,  it  will  be  used 
as  the  signal  source  to  test  the  implementation  of  the 
matched  filter.  The  matched  filter  requires  the  loading  of  a 


reference  code  sequence.  A  reference  code  generator  will  be 
Implemented  to  perform  this  function  and  is  basically  a  PN 
sequence  generator  and  timing  control  elements.  The  envelope 
detector,  accumulator,  and  decision  circuit  at  the  output  of 
the  matched  filter  will  be  modeled  separately.  At  this 
point,  the  modules  will  be  Integrated  for  implementation  of  a 
single  user  system.  The  additional  transmitters  necessary  to 
construct  a  multiple'-access  system  will  then  be  modeled. 

Pseudonoise  Sequence  Generator.  The  code  sequence 
produced  by  the  generator  depends  on  the  number  of  stages  of 
linear  shift  registers,  the  feedback  connections,  and  the 
initial  conditions.  For  this  system,  the  code  sequence  will 
be  generated  by  an  Hostage  linear  simple  shift  register 
(SSRG)  with  appropriate  feedback  to  produce  a  maximal  length 
sequence.  To  produce  the  maximal  length  sequence,  the 
feedback  connections  are  determined  by  the  selection  of  a 
primitive  11th  degree  polynomial  over  a  Galois  field  of  two. 
The  polynomial  selected, 

p(x)  «  x^^  +  x^  +  1  (41) 


will  produce  a  pseudonoise  sequence  with  a  repetition  period 
of  2^^  -  1  >  2047  code  symbols  (chips) .  Code  sequences 
with  long  repetition  periods  are  typically  selected  as  they 
result  in  code  sequences  that  are  virtually  decorrelated  for 
a  shift  of  a  single  chip.  A  high  degree  of  correlation 
between  an  undesired  code  and  a  receiver  reference  will  cause 
an  increase  in  the  receiver's  false  alarm  rate  and  possibly. 


false  receiver  synchronization.  Hie  block  diagram 
reprassntation  of  the  PN  sequence  generator  to  be  modeled  Is 
shotm  belov  in  Figure  24. 


For  efficiency  in  terms  of  simulation,  the  PN  sequence 
generator  will  be  implemented  with  a  reset  capability  for 
consistent  and  accurate  initiation  of  the  PN  sequence.  This 
will  eliminate  the  timing  uncertainty  of  initiating  the 
generator  by  closing  the  feedback  loop  and  loading  initial 
conditions  as  described  by  Smirnov  (3:68). 

The  fundamental  cell  used  in  the  construction  of  the 
generator  is  the  SSRG  Cell  with  a  reset  capability  as  shown 
in  Figure  25.  The  cell  provides  the  signal  and  the  running 
sum  inputs  and  outputs  necessary  for  the  multistage 
implementation.  The  initial  condition  for  each  stage  is  set 
at  the  Reset  Unit  Delay  module  and  the  feedback  connection  is 
established  by  setting  a  logical  true  condition  at  the 
Constant  Generator.  The  reset  signal  is  propagated  through 
the  logical  AND  gate  for  a  simultaneous  reset  of  the  cascaded 
cells.  The  complete  implementation  of  the  Resetable  PN 
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As  implemented  at  this  point,  the  generator  module  uses 
the  simulation  clock  and  only  generates  one  sequence  value 
per  sample  interval.  In  order  to  increase  flexibility  at 
simulation  time,  an  impulse  train  is  used  to  drive  the  chip 
rate  rather  than  the  simulation  clock. 

QPSK  Direct  Sequence  Transmitter.  The  quadriphase 
modulation  scheme  for  this  implementation  is  straightforward 
as  sho%m  in  Figure  27.  The  single  Resetable  PN  Generator  is 
used  to  modulate  both  the  in-phase  and  quadrature  channels. 
The  logical  true  and  false  outputs  of  the  generator  are 
converted  to  numerical  signal  levels  of  +1  and  -1 
respectively.  Operation  of  the  PN  generator  is  controlled  by 
the  enable  signal  from  the  Impulse  Train  module.  The 
parameter  that  establishes  the  rate  of  the  impulse  train  is 
exported  for  specification  at  simulation  time.  This 
parameter  is  termed  CHIP  RATE  and  establishes  chip  rate  of 
the  PN  generator  in  units  of  Hertz.  The  logical  false  input 
to  the  PN  generator  disables  the  reset  feature  of  the 
generator. 

The  data  sources  for  modulation  are  independent  random 
sources  into  each  channel.  The  numerical  signal  levels  are 
-fl  and  -1  and  the  signal  rate  is  identical  for  each  source. 
The  signal  rate  is  controlled  by  the  single  exported 
parameter  DATA  RATE.  The  random  sources  are  initialized  with 
different  seeds. 

The  channel  outputs  are  combined  into  a  complex  number 
and  modulated  on  a  complex  carrier.  The  frequency  produced 
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by  the  Complex  Spectral  Shifter  module  is  controlled  by  the 
>  exported  parameter,  CARRIER  FREQUENCY,  which  sets  the 

transmitter  output  frecjuency  in  units  of  Hertz. 

Matched  Filter  Channel .  As  shown  in  Figure  28,  the 
I  matched  filter  channel  consists  of  a  tapped-delay  line 

matched  filter,  envelope  detector,  accumulator  device,  and  a 
threshold  detector. 

'  Matched  Filter  Module.  Conceptually,  the 

implementation  of  the  matched  filter  is  in  the  form  of  a 
tapped  delay  line  followed  by  a  passive  filter  matched  to  a 
single  PN  chip  waveform.  Figure  29  illustrates  the  delay 
line  matched  filter  configuration  for  this  channel.  The 
delay  line  matched  filter  is  configured  to  recognize  a 
particular  code  sequence  of  finite  length  and  that  sequence 
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Figure  28.  Matched  Filter  Channel 
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Figure  29.  Tapped  Delay  Line  Matched  Filter  (11:70) 


only.  Each  delay  element  has  a  delay  equal  to  the  period  of 
the  expected  code  clock  so  that  each  element  contains  energy 
corresponding  to  only  one  code  chip  at  any  one  time  (13:30). 

The  matched  filter  is  a  passive  device  that  maximizes 
the  output  peak  signal-to-mean  noise  ratio.  The  signal  to  be 
detected  is  a  segment  of  the  input  signal  r(t)  of  length  T^ 
that  is  a  function  of  the  PN  code  modulation  only.  The 
impulse  response  of  the  filter  is  a  time-reversed  form  of  the 


signal  as  given  by  (11:69): 


h(t) 


r(Vt) 

0 


0  ^  t  £ 

o 

otherwise 


(42) 


or  in  terms  of  Fourier  transforms, 

H(w)  -  R*(w)  •e"^'^°  (43) 

Expressing  the  segment  length  as  ^  ^  encoded  PSK 

signal  takes  the  form, 

M 

r(t)  -  E  d  •p(t-(n-l)T  ]  (44) 

n=l  " 

where  «  {-1,+1)  is  the  polarity  of  the  n-th  chip  and 
p(t)  is  the  chip  pulse  shape  for  a  bandpass  implementation, 


P(t) 


yjp • cosw^t 
0 


0  s  t  <  T„ 
c 

otherwise 


(45) 


Taking  the  Fourier  transform  of  r(t)  and  substituting  its 
complex  conjugate  into  the  equation  for  H(w)  results  in  the 
implementation  of  Figure  29  and  the  following  expression 
(11:71) : 

H(w)  -  P*(w)  "  d  e-Jw(M-n+l)Tc  ,, 

n»l  " 
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Th«  simulation  model  takes  a  modified  form  of  the 
matched  filter  described  above.  Taking  the  inverse  Fourier 
transform  of  the  matched  filter  impulse  response  gives, 

»*  * 

h(t)  -  r  dp  t(M-n+l)T  -  t]  (47) 

n-1  "  ° 

Recognizing  that  d^  is  real,  an  equivalent  form  for  h(t)  is, 

M 

h(t)  -  r 
n=l 

The  simulation  model  establishes  the  tap  values  as  the 
complex  conjugate  of  d^,  given  below  and  eliminates  the 
filter  shown  in  Figure  29. 

dn.  -  d^  PC(M-n+l)T^,  -  tj  (49) 

The  modules  to  implement  this  filter  are  provided  by 
BOSS.  A  single  stage  of  the  matched  filter  is  shown  in 
Figure  30.  The  reference  sequence  is  loaded  into  the  top 
delay  channel  of  the  single  filter  cell.  A  load  enable 
signal  controls  the  sample  and  hold  module  and  locks  in  the 
reference  code  for  that  particular  matched  filter  cell  at  the 
end  of  the  load  sequence.  The  complex  conjugate  of  the 
reference  code  is  multiplied  with  the  incoming  signal  of  the 
bottom  delay  channel.  The  output  of  the  multiplier  is  added 
to  the  running  sum  output  of  the  set  of  cascaded  cells  that 
comprise  the  complete  matched  filter. 


dj^  P[(M-n+l)T^  -  t)] 
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Figure  30.  Correlate  Single  Stage 


The  complete  matched  filter  is  implemented  using  the 
Correlate  Multi-Stage  module  supplied  by  BOSS  as  shown  in 
Figure  31.  The  dashed  lines  indicate  use  of  the  replicate 
function.  The  replicate  function  allows  user-specification 
of  the  number  of  individual  stages.  The  length  of  the  filter 
is  set  at  simulation  time  with  the  parameter  I  OF  STAGES. 

The  maximum  length  of  the  filter  is  2047  stages.  The  initial 
running  sum  input  is  initialized  to  zero  and  the  delay  line 
outputs  are  terminated  with  sink  modules. 
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Envelope  Detector.  As  described  earlier,  the 
envelope  detector  consists  of  a  square-law  device  followed  by 
a  low-phss  filter.  A  simpler  Implementation  for  the 
simulation  model  will  provide  an  equivalent  output. 

The  sinusoidal  output  of  the  matched  filter  is 
amplitude-modulated  by  the  correlation  of  the  reference 
sequence  and  the  incoming  signal.  As  expressed  in  Eq  (26), 
the  matched  filter  output  is 

U(t)  -  P  MT^  Rpj,(t-MT^)  [dj(t)C0SW^t  +  dQ(t)8inw^t]  (26) 

The  desired  output  of  the  envelope  detector  given  in  Eq  (28) 
can  also  be  achieved  by  multiplying  the  incoming  signal  by 
its  complex  conjugate  and  reducing  the  gain  by  half.  This 
will  be  the  method  of  simulation  as  justified  in  the 
following  development. 

The  Incoming  signal  multiplied  by  its  complex  conjugate 

gives 

u(t)  u*(t)  «  (P  MT^  Rpjj(t-MT^)  |dj(t)cosw^t  +  dQ(t)sinw^tl^ 

where  |dj(t)cosw^t  +  dQ(t)8inw^t|  *  J2  .  With  the  bit 
energy  as  unity,  it  follows  that 

u(t)u*(t)  -  2[MRp^(t-MT^)]2  (50) 

This  result  is  twice  the  value  of  the  expected  signal  derived 
in  Eq  (28) .  It  is  therefore  valid  to  model  the  envelope 
detector  as  shown  in  Figure  32. 
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Figure  32.  Envelope  Detector 


Accumulator.  The  accumulator  circuit  performs  a 
running  sum  of  the  Incoming  signal  over  the  period  of  the  PN 
code.  As  shown  In  Figure  33,  a  cascade  of  delay  elements 
produce  the  necessary  delay.  The  size  of  both  multi-stage 
delays  Is  specified  by  the  parameter  0. 25*SAMPLES/PN  PERIOD. 


Figure  33.  Accumulator 


Threshold  Detector.  The  threshold  detector  must 
continuously  compare  the  Input  signal  to  a  specified 
threshold.  When  the  threshold  Is  exceeded,  the  detector  must 
transition  to  a  logical  true  condition  In  order  to  Initiate 
the  reset  function  In  the  PN  Sequence  Generator.  As  shown  in 
Figure  34,  the  comparison  value  Is  produced  by  the  Constant 
Generator  module  with  the  parameter  THRESHOLD.  When  the 
threshold  value  Is  exceeded,  the  comparator  circuit 
transitions  to  a  logical  true  condition  for  one  sample 
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Figure  34.  Threshold  Detector 


period.  The  Hold  True  Input  nodule  maintains  the  logical 
true  condition  through  the  remainder  of  the  simulation  as 
required  by  the  PN  Sequence  Generator. 

Quadrature  Correlator.  The  implementation  of  this 
simulation  model  follows  directly  from  the  block  diagram 
model  presented  in  Figure  7.  The  inputs  to  the  simulation 
model  shown  in  Figure  35  consist  of  the  received  signal  and 
the  local  PN  code-modulated  carrier. 
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Th«  real  and  imaginary  components  of  the  received  signal 
are  individually  multiplied  by  the  real  and  imaginary  outputs 
of  the  local  PN  code-modulated  carrier.  The  local  PN  code¬ 
modulated  carrier  is  a  complex  signal  with  the  real  and 
imaginary  components  corresponding  to  the  required  replica 
signals  for  the  in-phase  and  quadrature  channels  as  given 
below: 


Sj(t)  -  y2PPNj(t+r)cos(w^t+^)  (51) 

SQ(t)  -  y2P  PNQ(t+r)ain(w^t+^)  (52) 

The  integrate-and-dump  circuit  is  a  modified  version  of 
the  Integrate  And_Dump  module  provided  by  BOSS.  The 
modification  shown  in  Figure  36  provides  the  necessary  gain 
and  the  control  elements  to  output  the  integration  value  at 
the  end  of  a  every  PN  code  period.  The  output  of  the 
integrator  circuits  is  given  by  Eg  (16)  and  repeated  below: 


nT^ 

Xj(t)dt  -  Pdj(t) 
0 


nT_ 

c 

PN  j  ( t )  PN  j  ( t+T )  cos<p  [  l+cos2w^t  ]  dt 
0 


+  PdQ(t) 


.nT^ 

PNQ(t) PNj (t+T) sin^[ l-cos2w^t]dt 
0 


(16) 


The  gain  of  integrator  circuits  is  the  power  given  by  P. 

With  the  bit  energy  as  unity,  this  gain  is  equal  to  1/Tc. 

This  gain  is  provided  by  the  Gain  module  and  set  at 
simulation  time  by  the  CHIP  RATE  parameter  specified  also  for 
the  PN  Sequence  Generator.  The  Impulse  Train  module  selects 
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the  integrator  output  at  the  end  of  the  PN  code  period  and 
also  resets  the  integrator  for  the  start  of  the  next  period. 
At  all  other  times,  the  output  of  the  module  is  zero. 

Additicnal  Requirements.  Before  the  subsystem 

modules  can  be  integrated,  some  additional  modules  are 
necessary.  A  module  to  generate  the  seguence  for  the 
reference  channel  of  the  matched  filter  is  reguired  as  well 
as  a  module  to  generate  the  local  reference  for  the 
quadrature  correlator. 

Matghgd  Filter  Reference  Generator.  This  module 
loads  a  subsequence  of  the  pseudonoise  code  into  the 
reference  channel  of  the  matched  filter.  The  strategy  is  to 
use  the  last  M  code  chips  as  the  reference  sequence  where  M 
is  the  length  of  the  filter.  There  are  two  benefits  to  this 
method:  1)  the  matched  filter  generator  is  enabled  for  one 
full  code  period  and  the  appropriate  values  are  loaded 
regardless  of  the  filter  length  selected,  and  2)  by  using  the 
last  M  code  chips,  the  correlation  peak  occurs  at  the  instant 
of  the  start  of  a  new  PN  code  period.  This  prevents 
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initializing  the  local  PN  code  generator  with  different 
initial  conditions  for  different  Matched  filter  lengths.  The 
correlation  peak  always  signals  the  start  of  a  new  period. 

The  reference  generator  must  produce  a  replica  of  the  PN 
code-modulated  carrier.  To  accomplish  this,  the  locally 
generated  reference  PN  sequence  will  modulate  a  carrier  with 
the  same  timing  epoch  as  accomplished  in  the  transmitter. 

This  method  is  necessary  in  order  to  produce  an  exact 
reference  sequence  in  the  alternate  form  taken  for  the 
matched  filter. 

The  matched  filter  generator  is  shown  in  Figure  37.  The 
load  sequence  takes  place  in  the  first  PN  code  period  of  the 
simulation.  The  PN  generator  is  driven  by  the  impulse  train 
at  the  specified  chip  rate  and  then  modulates  a  carrier  at 
the  specif ed  frequency.  At  the  end  of  the  code  period,  the 
reference  generator  output  is  disabled.  The  reference  signal 
is  loaded  in  the  matched  filter  in  a  time-reversed 
orientation  and  locked  in  the  reference  channel. 

The  reference  generator  has  two  output  signals.  The 
REFERENCE  OUTPUT  signal  is  the  desired  reference  sequence. 

The  REFERENCE  LOAD  ENABLE  signal  is  a  logical  output  that 
enables  the  sample-and-hold  modules  of  the  matched  filter. 

The  timing  of  the  enable  signal  is  controlled  by  the  RESET 
SAMPLE  i  parameter  and  is  set  for  the  end  of  the  code 
sequence.  The  CHIP  RATE  and  CARRIER  FREQUENCY  parameters  are 
set  to  match  the  parameters  of  the  transmitter. 
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Figure  37.  Matched  Filter  Reference  Generator 


Lssal  Bfefier.gDSS  Generator.  This  module  produces 
the  local  reference  required  at  the  input  to  the  quadrature 
correlator.  This  module  is  similar  to  the  Matched  Filter 
Reference  Generator  described  above.  The  reference  signal 
for  the  quadrature  correlator  must  be  a  complete  complex 
signal.  The  real  and  imaginary  parts  provide  reference 
signals  for  the  in-phase  and  quadrature  channels  of  the 
correlator.  The  real  part  provides  the  PN  code-modulated 
cosine  signal  vhlle  the  imaginary  part  provides  the  modulated 
sine  function.  This  requirement  necessitates  the  structure 
shown  in  Figure  38.  The  PN  code  generator  feeds  both  sides 
of  the  Make  Complex  module.  In  the  matched  filter  case,  only 
the  real  signal  is  required  and  the  PN  code  generator  must 
only  feed  the  real  side  of  the  Make  Complex  module  (see 
Figure  37) .  The  reset  signal  to  the  Local  Reference 
Generator  is  used  to  produce  a  reset  strobe  to  the  PN 


Figure  38.  Local  Reference  Generator 


Sequence  Generator  and  select  the  PN  code>modulated  carrier 
for  output.  Once  again,  the  CHIP  RATE  and  CARRIER  FREQUENCY 
parameters  match  the  parameters  of  the  transmitter. 

Sinqlg-yggr  gY-stga  KQdsl 

The  single-'user  system,  shown  below  in  Figure  39,  is  an 
integration  of  the  subsystem  models  presented  so  far.  The 
single-user  system  is  discussed  in  further  detail  during  the 
simulation  and  analysis  of  the  next  chapter. 

Multiple-Access  System  Model 

In  code-division  multiple-access  (CDMA)  all  transmitters 
operate  at  the  same  nominal  frequency  and  simultaneously  use 
the  entire  bandwidth.  The  important  feature  of  CDMA  is  the 
absence  of  time  coordination  needed  between  the  various  users 
in  the  system.  Each  user  is  assigned  a  unique  PN  code  as  its 
spreading  function.  Assuming  synchronization  at  the 
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receiver,  the  intended  user's  spreading  function  will 
correlate  with  the  identical  despreading  function  of  the 
receiver.  Following  this  correlation,  only  the  original, 
modulated  waveform  remains.  Other  transmitted  waveforms 
appearing  at  this  receiver  will  not  be  correlated  and  are 
effectively  spread  and  appear  as  noise  to  the  demodulator 
(17:268-271) . 

The  correlation  noise  is  minimized  with  the  utilization 
of  sets  of  PN  sequences  with  low  cross-correlation 
properties.  As  seen  in  the  analytic  evaluation,  the 
reduction  of  the  cross-correlation  noise  power  and  associated 
increase  in  the  output  signal-to-noise  ratio  of  the 
correlation  receiver  is  an  important  factor  in  maximizing  the 
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total  number  of  simultaneous  users.  For  maximal  length 
sequences,  as  used  in  this  system,  only  small  sets  (five  or 
less)  can  have  good  cross-correlation  properties.  Large  sets 
of  maximal  length  sequences  generally  have  poor  cross¬ 
correlation  properties  and  are  typically  not  suitable  for  a 
large  number  of  users  in  a  CDMA  system  (17:281). 

For  the  multiple-access  simulation  model,  the  number  of 
other  users  is  set  at  four.  The  QPSK  Direct  Sequence 
Transmitter  will  act  as  the  target  transmitter.  The 
spreading  code  for  the  target  transmitter  is  the  maximal 
length  sequence  generated  by  the  PN  Sequence  Generator.  The 
transmitter  model  for  the  other  users  is  shown  in  Figure  40. 


Figure  40.  Multiple-Access  Transmitter 


This  transmitter  is  identical  to  the  QPSK  Direct 


Sequence  Transmitter  shown  in  Figure  27  with  the  exception  of 


the  source  for  the  spreading  code.  For  the  Multiple-Access 
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Transmitter,  the  Random  Data  module  provides  a  random  data 
stream  as  the  spreading  function.  It  has  been  shown  that 
purely  random  sequences  provide  slightly  better  CDMA 
performance  than  maximal  length  PN  codes.  This  difference 
can  be  made  negligibly  small  by  a  judicious  selection  of  the 
PN  sequences  used  (18:33.3.1).  Since  the  simulation  effort 
is  to  distinguish  the  target  transmitter  only,  use  of  the 
random  data  spreading  codes  is  further  justified. 

The  transmitter  set  for  the  multiple-access  system  is 
shown  in  Figure  41.  Carrier  frequencies  for  all  four 
transmitters  are  set  to  the  same  frequency  by  the  single 
parameter  TARGET  CARRIER  FREQUENCY.  The  seeds  for  the  random 
spreading  codes  of  each  Multiple-Access  Transmitter  are 


Figure  41.  Multiple-Access  Transmitter  Set 
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different  and  preset.  The  DATA  RATE  and  CHIP  RATE  parameters 
for  each  Multiple-Access  Transmitter  are  available  for 
specification  as  are  the  TARGET  DATA  RATE  and  TARGET  CHIP 
RATE  parameters  for  the  target  transmitter.  Additive  white 
Gaussian  noise  (AWGN)  is  specified  with  two  parameters.  The 
one-sided  bandwidth  is  set  with  the  parameter  NOISE  BW  and 
the  output  signal-to-noise  ratio  is  set  with  the  parameter 
SNR  (DB) .  The  Random  Start  modules  ensure  asynchronous 
operation  of  the  transmitters. 

The  multiple-access  system  is  shown  in  Figure  42. 
Parameter  designations  in  the  individual  modules  are  assigned 
to  facilitate  specification  at  simulation  time.  The  single 
specification  of  the  parameters  for  the  target  transmitter's 
carrier  frequency  and  chip  rate  also  specifies  associated 
parameters  in  other  modules.  The  parameter  TARGET  CARRIER 
FREQUENCY  sets  the  carrier  frequencies  for  the  other 
multiple-access  transmitters  and  sets  the  receiver's  local 
carrier  frequency  in  the  matched  filter  and  local  reference 
generators.  The  TARGET  CHIP  RATE  parameter  also  sets  the 
chip  rate  in  matched  filter  and  local  reference  generators  as 
well  as  the  gain  for  the  quadrature  correlator.  Simultaneous 
specification  of  these  parameters  reduces  the  simulation 
effort  and  avoids  possible  input  errors. 

Each  of  the  primary  modules-the  PN  Sequence  Generator, 
the  QPSK  Direct  Sequence  Transmitter,  the  Matched  Filter 
Channel,  and  the  Quadrature  Correlator  were  simulated  to 
confirm  operability.  These  simulations  also  indirectly 
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Figure  42.  Multiple-Access  System 

confirmed  operability  for  the  Matched  Filter  Reference 
Generator  and  Local  Reference  Generator  modules.  The 
simulation  results  are  documented  in  the  Appendix. 

In  the  next  chapter,  system  operability  will  be 
confirmed  through  simulation  of  the  Single-User  System. 
Further  simulations  will  then  examine  the  performance 
characteristics  of  interest. 
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£jl  Simulation 

Each  of  the  subsystem  simulation  models  described  in 
Chapter  IV  have  been  validated  with  methods  and  results 
documented  in  the  Appendix.  The  primary  purpose  of  the 
validation  is  to  verify  the  simulation  model  prior  to 
commitment  to  time-consuming  Monte-Carlo  simulations.  The 
length  of  the  validation  simulations  are  set  to  adequately 
verify  the  expected  end-to-end  behavior  of  a  system  or 
characterize  the  input/ output  behavior  of  the  subsystem 
modules. 

The  validation  mode  of  simulation  is  also  a  highly 
flexible  and  useful  tool  in  conceptual  communication  system 
studies.  It  develops  a  useful  and  realistic  appreciation  of 
actual  system  behavior  and  serves  as  a  useful  adjunct  to 
detailed  analysis  (8:58).  Development  of  a  model  that  will 
allow  a  realistic  appreciation  of  actual  system  behavior 
under  a  variety  of  scenarios  is  the  goal  of  the  simulation 
effort.  However,  a  simulation  model  must  be  proven  against 
theoretical  expectations  before  it  can  be  freely  exercised. 
Extensive  simulation  against  the  analytic  expectations 
developed  in  Chapter  III  is  the  immediate  concern.  Before 
proceeding  into  extensive  simulation,  it  is  necessary  to 
examine  the  system,  review  the  analytic  expectations,  and 
determine  which  system  performance  characteristics  can  be 
effectively  explored  via  simulation. 
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SimulatiOD  Approach 

The  mean  acquisition  time  is  the  performance 
characteristic  of  interest.  As  shown  in  Chapter  III,  the 
performance  of  the  matched  filter  as  a  function  of  input 
signal-to-noise  ratios,  matched  filter  length,  and  the  number 
of  other  users  in  the  system  directly  influences  the 
acquisition  time  performance.  Specifically,  the  matched 
filter  output  signal  must  be  tested  against  the  single  and 
multiple-pulse  detection  requirements  for  a  specified 
If  the  matched  filter  output  meets  the  specified  false  alarm 
rate,  the  incoming  PN  code  will  be  successfully  acquired  in  a 
known  time  frame.  The  accumulator  threshold  is  set  to 
process  the  required  number  of  correlation  peaks  out  of  the 
matched  filter  prior  to  initiating  the  local  PN  generator. 

The  number  of  correlation  peaks  necessary,  and  therefore  the 
acquisition  time,  is  specified  by  the  multiple-pulse 
detection  requirements  given  in  Figure  10. 

The  first  step  is  to  confirm  system  operability.  The 
single-user  input  is  used  to  validate  the  system  simulation 
model.  Following  validation,  a  series  of  simulations  will 
characterize  matched  filter  performance  for  both  the  single 
user  and  multiple-access  inputs.  Within  this  series  of 
simulations,  the  first  set  will  quantify  the  level  of  partial 
autocorrelation  noise  produced  by  the  matched  filter. 

The  next  set  of  simulations  will  determine  the  matched 
filter  peak  output  SNR  over  a  selected  range  of  input 
signal-to-noise  ratios  and  matched  filter  lengths.  From 
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these  results,  appropriate  threshold  levels  will  be 
determined  for  the  theoretical  single-pulse  or  multiple- 
pulse  detection  requirements.  Simulations  will  then  test  the 
matched  filter  against  the  detection  requirements  for  a 
specified  probability  of  false  alarm.  The  objective  is  to 
actively  test  matched  filter  performance  against  theoretical 
expectations.  If  successful,  the  simulation  model  is  proven 
suitable  to  simulate  realistic  operating  scenarios. 

System  Validation 

The  single-user  system  serves  as  the  base  simulation 
platform.  The  purpose  of  the  validation  with  the  single-user 
system  is  to  demonstrate  that  the  subsystem  models  will 
perform  as  expected  in  an  integrated  system.  The  single-user 
system  is  presented  below  in  Figure  43. 

A  single  transmitter  with  an  additive  white  Gaussian 
noise  (AWGN)  channel  is  implemented.  At  the  end  of  the 
initial  load  sequence,  the  Matched  Filter  Reference  Generator 
sends  a  ready  signal  to  select  the  transmitter  output.  The 
output  of  the  Quadrature  Correlator  is  terminated  into  a  SINK 
module-the  focus  of  the  simulation  effort  is  in  examining  the 
initial  acquisition  characteristics. 

Simulation  Parameters .  In  order  to  offer  flexibility  at 
simulation  time,  several  key  parameters  are  available  for 
specification.  For  the  transmitter,  the  parameters  DATA 
RATE,  CHIP  RATE,  and  CARRIER  FREQUENCY  are  set  in  units  of 
Hertz.  Note  that  the  one-time  specification  of  the  CHIP  RATE 
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and  CARRIER  FREQUENCY  parameters  not  only  sets  the  rates  for 
the  transmitter,  but  also  sets  rates  in  the  Matched  Filter 
and  Local  Reference  Generators.  Additionally,  the  CHIP  RATE 
parameter  establishes  the  gain  of  the  Quadrature  Correlator. 
The  noise  channel  has  parameters  to  specify  the  signal-to- 
noise  ratio  (SNR) ,  noise  bandwidth  (NOISE  BW) ,  and  the  time 
noise  is  activated  (NOISE  ON-TIME) . 

At  the  receiver,  the  matched  filter  is  specified  with 
the  parameter  #  OF  STAGES.  The  maximum  length  of  the  matched 
filter  is  2047  stages.  The  remaining  parameter  for  the 
Matched  Filter  Reference  Generator  specifies  the  number  of 
samples  in  one  PN  code  period.  This  parameter,  tSAMPLES/PN 
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PERIOD,  stops  the  load  sec[uence  of  the  matched  filter 
reference  channel,  locks  in  the  reference  values,  and  selects 
the  transmitter  output. 

The  accumulator  has  four  delay  modules  to  provide  a 
total  delay  equal  to  a  full  PN  code  period.  The  required 
input  for  the  parameter  0.25*SAMPLES/PERIOD  is  one-quarter  of 
the  number  of  samples  in  a  single  PN  period.  The  output  of 
the  accumulator  is  compared  against  a  threshold  value.  The 
threshold  is  set  with  the  parameter  ACCUMULATOR  THRESHOLD. 

Many  of  the  above  parameters,  such  as  RESET  SAMPLE  #  and 
NOISE  BW,  will  not  change  through  a  series  of  simulations. 
These  parameters  are  preset.  The  remaining  parameters  are 
exported  for  specification  at  each  simulation.  In  order  to 
confirm  the  system  model  operates  as  intended,  the  exported 
parameters  were  specified  as  shown  in  Table  2. 

Table  2.  Validation  Parameters 

STOP-TIME  ■  6.2E>3 

OT  «  6.25e-8 

SNR  (DB)  B  10 

NOISE  ON-TIME  a  6.2E-3 

ACCUMULATOR  THRESHOLD  m  60000.0 

A  OF  MF  STAGES  b  128 

DATA  RATE (HZ)  b  1200 

CHIP  RATE  (HZ)  b  2.0e6 

CARRIER  FREQUENCY  (HZ)  b  4 . 0e6 


Simulation  Results.  A  matched  filler  length  of  128 
stages  was  arbitrarily  chosen  for  the  validation  simulation. 
The  threshold  value  was  set  to  require  the  processing  of  four 
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correlation  peaks  (four  PN  code  periods)  at  the  output  of  the 
envelope  detector  output.  The  additive  white  Gaussian  noise 
(AW6N)  is  disabled  for  the  validation. 

During  the  first  PN  code  period  of  the  simulation,  the 
matched  filter  is  loaded  with  the  reference  sec[uence.  At  the 
end  of  the  next  period,  the  matched  filter  output  displays  a 
full  correlation  peak  as  sho%m  in  Figure  44.  With  the  length 
of  the  filter  equal  to  128,  the  amplitude  of  the  main 
correlation  peak  is  also  this  same  value. 


Figure  44.  Matched  Filter  Output 


The  output  of  the  accumulator  exceeded  the  threshold 
value  at  5.1175  milliseconds,  or  five  PN  code  periods  into 
the  simulation.  With  the  matched  filter  load  sequence 
occur ing  in  the  first  period,  the  accumulator  output  breaks 
the  threshold  value  of  60,000  in  the  expected  time  frame  as 
sho%m  in  Figure  45.  And  in  response,  the  output  of  the 
threshold  detector  moves  to  a  logical  true  value  as  shown  in 
Figure  46.  This  transition  initiates  the  Local  Reference 
Generator  as  shown  in  Figure  47. 
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Figure  45.  Accuaulator  Output 
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Figure  46.  Threshold  Detector  Output 
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Figure  47.  Local  Reference  Generator  Output 


The  Quadrature  Correlator  is  now  activated  for  full-period 
correlation  of  the  incoming  signals.  Figure  48  shows  the 


correlator  output  at  the  end  of  one  PN  code  period  following 
initiation  of  the  Local  Reference  Generator.  The  correlation 
peak  achieves  the  expected  maxinum  value  of  8.380418- 10^ 
indicating  synchronization  with  the  received  signal. 
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Figure  48.  Quadrature  Correlator  Output 

The  results  of  this  simulation  coupled  with  the 
subsystem  simulations  given  in  the  Appendix  validate  system 
operability  of  the  model.  The  next  task  is  to  confirm 
matched  filter  performance. 

Fi-Uer  Performance 

The  operating  characteristics  of  the  matched  filter  is 
the  key  to  system  performance.  For  a  given  signal-to-noise 
ratio  and  probability  of  detection,  the  matched  filter  output 
will  be  examined  to  determine  if  the  correlation  peaks  occur 
at  the  expected  value  for  the  specified  probability  of  false 
alarm  (P;^) . 

Determination  of  the  P^^  via  simulation  is  a  multi-step 
process.  First,  the  average  noise  power  at  the  output  of  the 
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matched  filter  is  determined  via  simulation  over  a  selected 
range  of  matched  filter  lengths  and  input  signal-to-noise 
ratios.  The  associated  peak  output  signal-to-noise  ratios 
are  then  determined  analytically.  The  peak  output  SNR 
distinguishes  whether  the  matched  filter  at  a  particular 
length  and  input  SNR  can  perform  at  single  or  multiple-pulse 
detection  requirements  for  a  specified  Appropriate 

threshold  values  are  then  set  to  confirm  matched  filter 
performance  at  the  specified 

Matched  Filter  Output  Noise  Power .  The  test  system  to 
determine  the  average  noise  power  for  the  single-user  system 
is  shown  in  Figure  49.  The  test  system  for  the  multiple 
access  system  uses  the  same  configuration  as  the  single-user 


Figure  49.  Matched  Filter  Noise  Test  System 


system  except  with  the  multi-transmitter  input.  The  output 
of  the  matched  filter  is  fed  to  the  Avg  Cmplx  Power  module. 
This  module  computes  the  average  output  power.  In  order  to 
get  an  accurate  measure  of  the  noise,  the  module  is  disabled 
while  a  main  correlation  peak  occurs  and  is  also  disabled 
during  the  initial  matched  filter  loading  period. 

The  results  of  the  simulations  are  given  in  Tables  3  and 
4  and  reflect  average  noise  power  given  an  input  signal  of 
unit  amplitude.  The  lengths  selected  for  the  matched  filter 
and  the  values  chosen  for  the  input  signal-to-noise  ratios 
will  be  consistent  throughout  subsequent  simulations. 


Table  3.  Single-User  Output  Noise  Power 


MF 

Length 

Target  Transmitter  SNR  Levels 

No  Noise 

3dB 

OdB 

-3dB 

-6dB 

64 

627.3 

731.4 

839.0 

1061.1 

1497.4 

128 

1305.1 

1515.7 

1733.6 

2163.4 

1 

3038.7 

256 

2610.9 

3034.6 

3478.6 

4377.4 

1 

6179.1 

512 

5231.3 

6081.7 

6945.1 

8675.0 

11922.3 

1024 

9658.9 

11280.4 

12952.9 

16314.7 

1 _ 

22723.5 

As  seen  in  the  first  column  of  Table  3,  the  partial 
autocorrelation  characteristics  of  the  PN  code  introduces  a 
significant  amount  of  noise  at  the  output  of  the  matched 
filter.  To  quantify  the  bias  introduced  by  this  noise,  the 
peak  output  SNR  must  be  determined.  The  peak  signal  power  is 
given  by  squaring  the  amplitude  of  the  main  correlation  peak. 


Table  4.  Multi'User  Output  Noise  Power 


Target  Transmitter  SNR  Levels 

No  Noise 

3dB 

OdB 

-3dB 

-6dB 

64 

2754.4 

2876.6 

2983.1 

3143.8 

3528.0 

128 

5816.7 

6043.3 

6252.0 

6651.5 

7454.6 

256 

11694.6 

12323.7 

12822.1 

13388.2 

15036.9 

512 

24024.8 

25130,7 

26042.1 

27076.8 

30145.2 

1024 

44973.1 

46304.1 

47731.9 

50623.3 

56452.8 

As  shown  in  Chapter  III,  the  expected  value  of  this  amplitude 
is  equal  to  the  length  of  the  filter.  This  method  provides 
the  peak  output  SNR  results  given  in  Tables  5  and  6.  Recall 
that  the  peak  output  SNR  is  the  key  element  in  determining 
acquisition  time  performance. 


Table  5.  Single-User  Peak  Output  SNR  (dB) 


MF 

Length 

Target  Transmitter  SNR  Levels 

No  Noise 

3dB 

OdB 

-3dB 

-6dB 

64 

8.15 

7.48 

6.89 

5.86 

4.37 

128 

10.99 

10.34 

9.75 

8.79 

7.32 

256 

13.99 

13.34 

12.75 

11.75 

10.26 

512 

17.00 

16.35 

15.77 

14.80 

13.42 

1024 

20.36 

19.68 

19.08 

18.08 

16.64 

The  partial  autocorrelation  of  the  entire  PN  code 
against  a  segment  of  itself  effectively  introduces  a  bias  to 
the  input  SNR.  In  other  words,  an  input  SNR  equal  to  3  dB 


85 


Table  6.  Multi-User  Peak  Output  SNR  (dB) 


MF 

Length 

Target  Transmitter  SNR  Levels 

No  Noise 

3dB 

OdB 

-3dB 

-6dB 

64 

1.72 

1.53 

1.38 

1.15 

0.65 

128 

4.50 

4.33 

4.18 

3.92 

3.42 

256 

7.48 

7.26 

7.09 

6.90 

6.39 

512 

10.38 

10.18 

10.03 

9.86 

9.39 

1024 

13.68 

13.55 

13.42 

13.16 

12.70 

appears  at  the  natched  filter  output  as  a  lower  input  SNR  due 
to  the  introduction  of  this  partial  autocorrelation  noise. 
Determination  of  this  bias  allows  accurate  comparison  of  the 
simulation  results  to  the  expected  performance.  The  bias  is 
determined  using  the  results  listed  in  the  first  column  of 
Table  5.  The  gain  of  the  matched  filter  is  equal  to  the 
length  of  the  filter  and  is  used  to  determine  the  effective 
input  SNR.  These  results  and  the  effective  input  SNR  due  to 
the  multi-user  influence  are  given  in  Table  7. 


Table  7.  Input  SNR  Bias 


MF 

Length 

Effective  Input  SNR  (dB) 

Correlation  Noise 
Only 

Correlation  Noise 
w/  Multi-User  Input 

64 

-9.91 

-16.34 

128 

-10.08 

-16.57 

256 

-10.09 

-16.60 

512 

-10.09 

-16.71 

1024 

-9.74 

-16.43 
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These  levels  establish  the  upper  bounds  of  input  SNR  levels 
for  which  the  siaulation  results  can  be  conpared  against  the 
analytic  expectations.  The  remaining  peak  output  SNR  levels 
given  in  Tables  5  and  6  are  used  to  distinguish  which  matched 
filter  lengths  and  associated  input  SNR  levels  can  perform  at 
the  single-pulse  detection  requirements.  The  minimum  peak 
output  SNR  levels  for  specified  probabilities  of  false  alarm 
were  given  in  Figure  8  and  summarized  in  Table  1  for  a 

»  0.9.  The  chosen  for  simulation  is  10  and  requires 
a  minimum  peak  output  SNR  of  13.7  dB  to  meet  single-pulse 
detection  requirements.  This  false  alarm  rate  can  be 
simulated  with  reasonable  time  and  memory  requirements. 

Hatgh.?d  ZUXsr  False  Alarm  Performance.  The  false  alarm 
probability  is  the  probability  that  a  threshold  crossing  will 
occur  during  the  interval  of  observation  when  the  receiver 
input  contains  only  noise.  If  a  false  alarm  occurs  in  the 
matched  filter  channel,  the  local  reference  generator  will  be 
improperly  initiated  and  the  quadrature  correlator  will  not 
indicate  synchronization.  The  penalty  for  this  false  alarm 
is  time-the  full  period  correlation  sequence  of  the 
quadrature  correlator  and  the  time  expended  by  the  matched 
filter.  If  the  matched  filter  is  in  a  multiple-pulse 
detection  sequence,  the  tine  becomes  significant. 

In  order  to  confirm  P^^  performance,  the  threshold  level 
must  be  accurately  determined.  There  are  two  methods  to 
confirm  performance  under  single-pulse  detection 
requirements.  The  first  method  sets  the  threshold  level  for 
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the  eignal  directly  at  the  output  of  the  Hatched  filter.  The 
governing  equation  for  this  method  is  (16:204): 


where  p  *  .  The  threshold  value  (a)  necessary  to  test 

for  a  particular  is  dependent  on  the  average  noise  power 
(o^  )  at  the  sampling  instant. 

The  second  method  establishes  the  threshold  level  at  the 
output  of  the  envelope  detector.  A  table  of  bias  levels 
6  »  f(m,p)  corresponding  to  the  specified  probability  of 
false  alarm  (10~^)  as  a  function  of  the  number  of  processed 
pulses  (m)  exist  for  1  6  p  6  12  and  1  6  m  S  150  for 
receivers  using  square-law  detectors  (19:38-45).  Given  the 
bias  level  from  the  table  and  ,  the  threshold  level  is 
determined  by: 


o  -  20^*6 


This  method  is  also  suitable  for  the  multiple-pulse  detection 
tests.  For  m  »  1  and  not  included  in  the  table  of  bias 
levels,  the  threshold  level  can  be  computed  using  the 
following  equation  (16:245): 


1  -  I 


2 


m-1 


where  I(x,y)  is  a  form  of  the  incomplete  Gamma  function. 
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Threshold  Settings.  The  first  aethod  will  be  used 
to  set  the  thresholds  for  the  single-pulse  detection  tests. 
This  method  is  sinple  and  direct  and  offers  a  certain  clarity 
in  terns  of  the  threshold  value  against  the  expected  peak 
signal  value.  The  second  method  was  the  only  method  found  to 
set  threshold  values  for  the  multiple-pulse  detection  tests. 

The  computed  threshold  settings  are  given  in  Tables  8 
and  9.  The  asterisk  (*)  indicates  a  threshold  setting  for 
single-pulse  detection  requirements.  If  a  peak  output  SNR 
falls  below  the  minimum  required  for  single-pulse  detection 
but  the  single-pulse  detection  threshold  is  less  than  the 
expected  peak  signal  amplitude,  a  single-pulse  detection  test 
will  be  run.  These  situations  are  indicated  by  a  double 
asterisk  (**).  For  the  multiple-pulse  detection  thresholds, 
the  niunber  of  pulses  necessary  is  indicated  in  parenthesis. 


Table  8.  Single-User  Threshold  Levels  “  lo”^) 


MF 

Length 

Target  Transmitter  SNl 

1  Levels 

No  Noise 

3dB 

OdB 

-3dB 

-6dB 

64 

18.56e3 

28.71e3 

32.94e3 

57.36e3 

113. 9e3 

(5) 

(8) 

(8) 

(13) 

(21) 

128 

24.11e3 

27.99e3 

38.93e3 

56.52e3 

109.77e3 

(2) 

(2) 

(3) 

(4) 

(7) 

256 

189.92 

204.75 

219.22 

245.92 

138.77e3 

* 

** 

** 

** 

(3) 

512 

268.84 

289.86 

309.76 

346.19 

405.85 

• 

* 

* 

* 

** 

1024 

365.30 

394.77 

423.02 

474.76 

560.30 

* 

# 

* 

* 

* 

Table  9.  Multi-User  Threshold  Levels  *  10~^) 


Multi-User  Input  SNR  Levels 

No  Noise 

3dB 

OdB 

-3dB 

-6dB 

64 

458.35e3 

(57) 

506. 3e3 
(61) 

539. 3e3 
(63) 

620. 7e3 
(70) 

845. 7e3 
(88) 

128 

426.96e3 

(20) 

475. 9e3 
(22) 

508. 9e3 
(23) 

593. 8e3 
(26) 

742. 5e3 
(30) 

256 

459.04e3 

(8) 

483. 7e3 
(8) 

503. 3e3 
(8) 

566. 5e3 
(8) 

725. 8e3 
(11) 

512 

443.18e3 

(2) 

564. 4e3 
(3) 

584. 9e3 
(3) 

608. Ie3 
(3) 

787. 5e3 
(4) 

1024 

788.24 

** 

799.82 

** 

812.10 

** 

836.29 

** 

883.12 

** 

The  next  step  is  to  determine  how  long  the  simulations 
must  run  in  order  to  confirm  performance  at  a  particular 
probability  of  false  alarm. 

fAlSfi  Alum  Performance  Results.  As  stated  in 
Chapter  II,  the  rule  of  thumb  is  to  run  the  simulation  long 
enough  to  count  at  least  ten  errors.  For  a  specific  P^^r  the 
period  of  simulation  is  equal  to  ^^/Pf^i  times  the  interval 
over  which  a  false  alarm  can  occur.  The  interval  over  which 
a  false  alarm  can  occur  is  different  for  single-pulse  and 
multiple-pulse  detection  sequences. 

For  the  single-pulse  detection  sequences,  a  false  alarm 
can  occur  over  any  two-chip  period.  The  matched  filter 
continuously  correlates  its  reference  sequence  with  the 
incoming  signal  and  the  main  correlation  peak  occurs  over  a 
two-chip  period.  During  all  other  two-chip  intervals,  a 
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correlation  peak  exceeding  the  specified  threshold  can  occur 
due  to  noise  or  the  cross*>correlatlon  characteristics  of  the 
Input  signal.  Therefore,  the  length  of  all  simulations 
testing  for  single-pulse  detections  Is  set  according  to: 


STOP  TIME 


DT- 


Samples 

Chip 


■  (2  Chips) 


(56) 


The  time  at  which  the  simulation  Is  stopped  Is  designated  by 
the  parameter  STOP  TIME  and  the  Interval  between  sampling 
Instants  Is  DT-the  discrete  time  Interval.  The  discrete  time 
Interval  Is  set  such  that  8  to  16  samples  occur  per  PN  code 
chip. 

The  Interval  over  which  a  false  alarm  can  occur  for 
multiple-pulse  detection  sequences  Is  the  number  of  pulses  to 
be  processed  times  an  entire  PN  code  period.  The  threshold 
value  for  a  multiple-pulse  detection  sequence  Is  crossed  by 
the  accumulated  amplitude  of  the  number  of  processed  pulses. 
The  length  of  all  simulations  testing  P^^  for  multiple-pulse 
detections  Is  set  according  to  the  equation: 


10 


STOP  TIME 


DT 


fa 


Samples 

Chip 

#  Chips 


PN  Period 


(«  Pulses)  (57) 


Obviously,  the  length  of  simulations  to  characterize 
multiple-pulse  detection  requirements  are  excessive. 


For  the  single-pulse  detection  sequences,  the  matched 
filter  channel  will  be  characterized  for  a  equal  to  10~^. 
Simulations  at  this  can  be  accomplished  in  5-6  hours  per 
table  entry.  Simulations  at  lover  values  of  P^^  become 
excessive  in  terms  of  time  and  memory  requirements.  A 
validation  simulation  will  be  run  for  the  multiple-pulse 
detection  sequences.  In  other  words,  the  simulation  will 
verify  that  the  specified  threshold  is  crossed  in  the 
expected  time  frame.  The  time  required  to  verify  multiple- 
pulse  detection  performance  at  a  reasonable  P^^  extends  well 
beyond  the  gain  in  terms  of  insight  to  system  performance. 

As  discussed  in  Chapter  II,  simulations  of  this  nature 
are  an  exercise  in  counting  errors.  The  test  system  for  the 
single-pulse  detection  requirements  is  shown  in  Figure  50. 
This  test  system  simply  counts  the  number  of  threshold 
crossings  over  the  simulation  period.  Two  threshold  crossing 
counts  are  taken.  The  first  count  tracks  the  number  of 
threshold  crossings  during  the  period  which  only  noise  is 
present.  This  count  provides  the  numbers  necessary  for  P^^ 
calculations.  The  second  count  provides  the  total  number  of 
crossings.  This  number  is  used  as  a  check  to  insure  the 
signal  is  detected  when  present. 

The  test  system  to  validate  the  multiple-pulse  detection 
requirements  is  shown  in  Figure  51.  In  this  system,  the 
threshold  value  is  set  according  to  values  in  Tables  8  and  9. 
When  the  threshold  is  crossed,  the  simulation  will 
automatically  terminate  and  the  time  at  which  the  threshold 
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Is  crossed  will  be  recorded.  This  time  period  Is  then 
Batched  against  the  expected  tiae  franes.  The  single-user 
input  depicted  in  each  of  the  above  test  systems  is 
substituted  with  the  Bulti-user  input  as  necessary. 

The  results  of  the  siaulations  are  given  in  Tables 
10  and  11.  In  the  single-pulse  detection  tests,  the  result 
is  a  value  for  false  alara  performance.  The  multiple-pulse 
detection  tests  are  validation  simulations  which  yield  pass 
or  fail  results. 

Table  10.  Single-User  Results 

MF  Target  Transmitter  SNR  Levels 

Length  - , - 1 - 1 - 1 - 


No  Noise 

3dB 

OdB 

-3dB 

-6dB 

64 

Pass 

Fail 

Fail 

Fail 

Fail 

128 

Pass 

Pass 

Pass 

Pass 

Pass 

256 

0.69e-3 

0.50e-3 

0.20e-3 

<1.0e-3 

Pass 

512 

1.19e-3 

1.09e-3 

0.89e-3 

0.69e-3 

0.30e-3 

1024 

1.39e-3 

1.39e-3 

1.39e-3 

1.09e-3 

0.89e-3 

Table 

LI.  Multi- 

-user  Pj^  I 

Results 

MF 

Length 

Target  Transmitter  SNR 

Levels 

No  Noise 

3dB 

OdB 

-3dB 

-6dB 

128 

Fail 

Fail 

Fail 

Fail 

Fail 

256 

Pass 

Pass 

Pass 

Fail 

Fail 

512 

Pass 

Pass 

Pass 

Pass 

Pass 

1024 

0.20e-3 

0.20e-3 

0.20e-3 

0.20e-3 

0.20e-3 

For  a  aatchad  filtar  of  langth  aqual  to  64  stages,  the 
■ultipla-pulse  detection  tests  failed  with  the  single-user 
input.  These  failures  are  attributed  to  the  high  level  of 
partial  autocorrelation  noise  produced  by  the  selected  code. 
These  correlation  peaks  add  linearly  froa  period  to  period. 
This  noise  coupled  with  the  AWGN  consistently  produced 
failure  indications  even  with  threshold  level  and  time  period 
adjustments.  This  filter  length  was  therefore  not  simulated 
with  the  multiple-user  input.  At  shorter  matched  filter 
lengths,  performance  becomes  increasingly  sensitive  to  the 
correlation  noise  produced  by  the  selected  PN  code.  It  is 
reasonable  to  ass\ime  that  a  PN  code  with  lover  partial 
autocorrelation  characteristics  would  produce  favorable 
results  at  these  matched  filter  lengths. 

In  the  multi-user  system,  code  selection  is  a 
significant  consideration.  Ideally,  a  set  of  orthogonal 
codes  is  used.  In  the  simulation,  the  effort  did  not  entail 
a  code  selection  process  for  the  multi-user  system-detection 
of  the  single  target  transmitter  was  the  objective. 

Therefore,  random  sequence  generators  substituted  controlled 
code  generation  by  the  other  users.  The  impact  was  increased 
cross-correlation  noise  and  higher  threshold  levels.  The 
effect  was  significant  at  a  matched  filter  of  length  equal  to 
128  stages.  Once  again,  it  is  reasonable  to  assume  that 
careful  code  selection  would  produce  favorable  results  at 
these  matched  filter  lengths. 
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Th«  remaining  results  confirm  the  expected  performance 
characteristics.  All  single-pulse  detection  tests  were 
successful.  The  validation  tests  for  multiple-pulse 
detections  were  also  successful  at  the  longer  matched  filter 
lengths.  Sensitivity  to  the  correlation  noise  at  these 
matched  filter  lengths  was  not  a  factor.  It  is  highly 
probable  that  the  successful  validation  simulations  would 
also  prove  successful  with  the  longer  simulations  necessary 
to  confirm  false  alarm  performance.  Recognizing  the  impact 
of  code  selection  and  the  correlation  noise  produced  by  the 
selected  PN  code,  the  simulation  model  is  suitable  to  provide 
a  realistic  appreciation  of  actual  system  behavior  under  a 
variety  of  scenarios. 
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The  prlnary  objective  of  this  study  has  been  to 
investigate  the  performance  characteristics  of  a 
synchronization  system  that  uses  an  auxiliary  matched  filter 
channsl  to  speed  acquisition  time.  A  quantitative  analysis 
was  conducted  to  characterize  mean  acquisition  time  under 
varying  conditions  of  matched  filter  length,  noise  levels, 
and  number  of  users  in  the  system.  A  simulation  model  of  the 
synchronization  system  was  constructed  using  the  Block 
Oriented  Systems  Simulator  (BOSS) .  The  objective  for 
simulation  was  to  confirm  the  analytic  expectations  and 
provide  a  proven  simulation  model  for  further  system 
analysis. 

The  quantitative  results  indicate  that  the  auxiliary 
matched  filter  channel  Improves  acquisition  time.  For 
signal-to-noise  ratios  above  OdB  and  a  particular  matched 
filter  length,  the  mean  acquisition  time  was  directly  a 
function  of  the  number  of  users  in  the  system.  Also,  mean 
acquisition  time  was  significantly  improved  for  signal-to- 
noise  ratios  above  OdB.  As  one  example,  for  a  matched  filter 
length  of  128  stages  and  64  users  in  the  system,  the  matched 
filter  channel  provides  a  gain  in  acquisition  time  of  156 
over  the  active  correlator  operating  alone.  As  the 
signal-to-noise  ratio  decreases  below  OdB,  the  variation  in 
mean  acquisition  time  caused  by  the  number  of  users  rapidly 
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declines.  For  signal-to-noise  ratios  below  -18dB,  mean 
acquisition  time  is  almost  directly  related  to  the  noise 
level-irrespective  of  the  number  of  users.  Mean  acquisition 
time  asymptotically  approaches  a  doubling  in  time  for  each 
1.5  dB  decrease  in  the  signal-to-noise  ratio. 

The  simulation  results  confirmed  the  analytic 
expectations  for  false  alarm  performance  of  the  matched 
filter  channel.  All  single-pulse  detection  requirements  at 
the  specified  were  favorably  confirmed.  Multiple-pulse 
detection  requirements  were  favorably  validated  with  shorter 
simulations.  The  false  alarm  interval  for  the  multiple-pulse 
detection  requirements  was  too  lengthy  to  confirm  performance 
at  the  specified 

Exceptions  to  the  multiple-pulse  detection  results  were 
noted  due  to  the  high  partial  autocorrelation  noise  produced 
by  the  selected  PN  code.  The  noise  generated  at  the  output 
of  the  matched  filter  by  this  correlation  of  the  selected  PN 
code  and  its  reference  segment  was  significant.  The  shorter 
matched  filter  lengths  proved  sensitive  to  the  linearly 
additive  qualities  of  the  correlation  noise  under  the 
multiple-pulse  detection  tests. 

The  simulated  probability  of  false  alarm  was  restricted 
to  10  .  Simulation  times  for  lower  false  alarm  rates  were 

excessive.  The  BOSS  provided  a  highly  flexible  simulation 
environment.  The  system  model  includes  flexibilities  to 
adjust  matched  filter  lengths,  noise  levels,  and  threshold 
levels.  With  a  minimum  of  effort,  a  different  pseudonoise 


code  may  be  selected,  the  modulation  scheme  adjusted,  and  the 
number  of  users  in  the  system  can  be  changed.  In  this 
effort,  the  simulation  model  was  used  to  statically  confirm 
the  analytic  expectations.  The  successful  static  tests  prove 
suitability  of  the  model  to  accurately  simulate  dynamic 
applications  of  system. 

Regonmgndationa 

The  following  recommendations  are  proposed  for  further 
study : 

1.  h  number  of  proposed  systems  use  a  similar  configuration 
to  speed  code  acquisition  (12;  14;  20;  21;  22).  None  of 
these  systems,  however,  implement  an  accumulator  circuit. 
Additionally,  Smirnov  does  not  describe  dynamic  operation 
of  the  accumulator.  Methods  to  control  the  accumulator 
should  be  explored.  Control  modules  for  the  accumulator 
can  be  added  to  the  present  simulation  model  to 
dynamically  simulate  and  analyze  system  performance. 

2.  A  standard  Monte  Carlo  technique  was  used  to  test  false 

alarm  performance.  For  probability  of  false  alarm  values 
below  10  ,  the  simulation  times  were  excessive.  A 

modified  Monte  Carlo  technique  such  as  importance 
seuapling  could  significantly  reduce  simulation  times  and 
allow  performance  evaluations  at  more  desirable  false 
alarm  rates.  Modified  Monte  Carlo  techniques  should  be 
Investigated  and  implemented  for  the  BOSS. 
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Appendix:  siattlatioD  Mgdsi  Validatigng 


The  purpose  of  this  appendix  is  to  validate  the 
simulation  models  presented  in  Chapter  III.  The  emphasis  in 
the  validation  phase  is  to  show  that  the  models  operate  as 
intended.  The  simulation  parameters  are  set  at  the  minimum 
necessary  to  confirm  expected  input/ output  behavior. 

Eg^VtdPnpjgg  Sequence  Generator 

The  first  test  of  the  PN  generator  confirms  the  length 
and  character  of  the  code.  The  reset  feature  of  the 
generator  is  then  tested.  The  configuration  for  both  tests 
is  shown  in  Figure  52. 


Figure  52.  PN  Generator  Test 


An  impulse  train  controls  the  chip  rate  of  the  PN 

generator.  The  system  parameter,  CHIP  RATE,  is  set 

throughout  the  validations  at  2  MHz.  The  simulation  is  set 

for  eight  samples  per  chip  period  which  requires  the  discrete 

_  —8 

time  interval  parameter,  DT,  to  be  set  at  6.25- 10  seconds 
per  sample.  For  the  first  test,  the  length  of  the  simulation 
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is  set  to  run  through  two  coaplete  code  periods  requiring  a 
STOP  TIME  equal  to  2.147- 10~^  seconds.  The  reset  feature  is 
not  of  interest  in  this  test.  The  occurence  of  the  reset 
strobe  is  set  to  the  end  of  the  simulation  period.  The 
simulation  parameters  are  given  in  Table  12. 


Table  12.  PN  Generator  Validation  Parameters 


STOP-TIME  ■  2.1470001E-3 
DT  a  6.25e-8 
RESET  SAMPLE  A  a  17970 
RESET  LENGTH  (ASAMPLES)  a  2 
CHIP  RATE  <HZ)  a  2.0e6 


The  length  of  the  code  is  2^^-l  ■  2047  chips  .  At  the  given 
chip  rate,  the  code  should  repeat  every  1.0235  milliseconds. 
The  signal  plots  of  Figures  53  and  54  were  taken  at  the  probe 
placed  at  the  input  to  the  sink  module. 


Figure  53.  PN  Generator  Output 


The  signal  plot  of  Figure  53  displays  the  code  at  the  start 
of  the  simulation.  The  dashed  presentation  of  the  signal 


101 


plot  was  produced  by  overlaying  the  plot  of  Figure  54  which 
starts  at  1.0235  Billiseconds.  The  dashed  line  confirms  that 
the  code  repeats  as  expected. 


Figure  54.  PN  Generator  Output 


The  autocorrelation  of  the  PM  generator  output  is  shown 
in  Figure  55  and  at  a  code  period  later  in  Figure  56.  As 
expected,  the  width  of  the  autocorrelation  waveform  is  twice 
the  chip  period  and  equal  to  one  microsecond.  The  second 
correlation  peak  occurs  at  exactly  one  code  period  after  the 
start  of  the  sequence  (1.0235  microseconds).  The  amplitude 
of  the  correlation  peaks  do  not  reach  unity  due  to  the 
sampling  rate  and  the  memory  limitations  of  this  particular 
post-processor  function.  The  correlation  function  will 
correlate  a  maximum  of  8000  samples.  The  plots  provided  use 
3435  samples  for  the  time  frame  to  be  correlated  and  3435 
samples  for  the  reference  sequence.  The  magnitude  of  the 
autocorrelation  is  the  length  of  the  correlation  Interval 
times  the  chip  period  which  is  equal  to  2. 146875  lo'^  as 
indicated. 
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The  reset  feature  is  now  tested  by  setting  the  reset 
strobe  to  occur  at  250  Billiseconds.  As  shown  in  Figure  57, 
the  reset  output  transitions  to  a  logical  true  condition  for 


two  sample  periods. 


Figure  57.  Reset  Strobe 


The  strobe  resets  the  PN  Generator  to  the  specified  set  of 
initial  conditions.  Figure  58  shows  the  initiation  of  a  new 
code  sequence  in  response  to  the  reset  signal. 


Figure  58.  PN  Generator  Output 


As  a  final  check,  the  magnitude  spectrum  of  the  PN  Generator 
output  is  given  in  Figure  59. 
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Power  (dB) 


Bvorogo  Power  1.0 
OdB  Ref.  Power  0.004145 


Figure  59.  Magnitude  Spectrun  of  the  PN  Output 


Q£S£  Pirjeg-t  Transmitter 

The  transmitter  output  vill  be  tested  to  confirm  the 
carrier  is  properly  modulated  by  the  data  and  PN  codes.  The 
test  system  shown  in  Figure  60  is  simply  the  transmitter 
with  its  output  terminated  into  a  sink  module.  Simulation 
probes  are  set  to  observe  the  transmitter  output  as  well  as 
the  data  and  PN  code  inputs. 

In  order  to  present  reasonable  detail  of  the  sinusoidal 
output,  the  sampling  frequency  is  set  for  sixteen  samples  per 
carrier  period.  The  carrier  frequency  is  set  at  4  MHz 
throughout  the  validations.  The  length  of  the  simulation  is 
influenced  by  the  data  rate.  The  STOP  TIME  is  set  to  observe 
at  least  one  data  bit.  Simulation  parameters  are  given  in 
Table  13. 


r 


5 


Figure  60.  Transmitter  Test  System 


Table  13.  Transmitter  Validation  Parameters 


STOP-TIME  ■  8.4999996E-4 
DT  s  1.5625e-6 

CARRIER  FREQUENCY  <HZ)  s  4 . 0e6 
CHIP  RATE  (HZ>  ■  2.0e6 
DATA  RATE (HZ)  ■  1200 


The  transmitter  output  shown  in  Figure  61.  The  plot  is 
overlayed  with  the  PN  code  over  the  same  interval  shown  in 
Figure  62.  Carrier  phase  reversals  occur  coincident  with  the 
the  PN  code  transitions.  Over  this  time  interval,  a  data  bit 
transition  also  occurs.  The  phase  reversal  due  to  this 
transition  is  sho«m  with  an  overlay  of  the  data  output  in 
Figure  63 . 


Figure  62.  PN  Code  Output 
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previously  in  Figure  39.  Probes  are  placed  at  the  outputs  of 


the  Batched  filter,  envelope  detector,  and  accumulator 
modules.  The  validation  simulations  were  run  without  noise 
and  set  to  run  through  three  full  PN  code  periods.  The  first 
period  is  dedicated  to  loading  the  matched  filter  reference 
sequence.  The  simulation  parameters  are  given  in  Table  14 
and  the  output  of  the  matched  filter  at  the  end  of  the  second 
PM  code  period  is  shovm  in  Figure  65. 


Table  14.  Matched  Filter  Validation  Parameters 


STOP-TIME  e  3.1E-3 

OT  e  6.25e-8 

SNR  (OB)  a  10 

NOISE  ON-TIME  s  3.1E-3 

ACCUMULATOR  TNRESHOLO  e  16000 

«  OF  MF  STAGES  s  128 

OATA  RATE (HZ)  s  1200 

CHIP  RATE  (HZ)  s  2 . 0e6 

CARRIER  FREQUENCY  (HZ)  s  4 . 0e6 


Figure  65.  Matched  Filter  Output  (Length  «  128) 


Figure  65  confirms  that  the  expected  amplitude  of  the  main 
correlation  peak  is  equal  to  the  length  of  the  filter.  The 
output  of  the  matched  filter  is  envelope-detected  through  a 
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square-law  device  and  a  low-pass  filter.  The  expected 
anplitude  of  the  sain  correlation  peak  is  the  square  of  the 
■atched  filter  length  or  (128)^  «  16,384  .  The  output  of 
the  envelope  detector  is  sho%m  in  Figure  66. 


Figure  66.  Envelope  Detector  Output 


The  output  of  the  Accumulator  is  sho%m  in  Figures  67 
and  68  for  the  end  of  the  second  and  third  PN  code  periods 
respectively.  Note  that  the  accumulated  amplitude  of  the 
sain  correlation  peak  in  Figure  68  is  twice  that  shown  in 
Figure  67.  Also  note  that  the  noise  due  to  the  partial 


Figure  67.  Accumulator  Output 
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X  10**3 


Figure  68.  Accuaulator  Output 


autocorrelation  of  the  selected  code  adds  linearly  from  one 
period  to  the  next-this  noise  is  not  random.  Similar  results 
were  obtained  for  matched  filter  lengths  equal  to  64,  256, 
512,  and  1024. 

Quadrature  Correlator 

The  Quadrature  Correlator  serves  a  specific  purpose  in 
the  simulation  effort.  Once  the  Matched  Filter  Channel  has 
initiated  the  Local  Reference  Generator,  it  is  the  correlator 
output  that  indicates  successful  synchronization  to  the 
incoming  signal.  The  amplitude  of  the  correlator  output 
signal  was  determined  in  Chapter  III  and  is  repeated  below: 


Z^{t) 


"  T  i 

nRpjj(T) 


(20) 


where  n  ■  2047  is  the  number  of  chips  in  a  complete  PN  code 
period.  When  synchronization  occurs,  Ro»(r}  -  1  and  Z„(t) 
reaches  its  maximum  value  of  8.380418* 10^. 


4 


4 
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The  purpose  of  the  validation  sinulation  for  the 
Quadrature  Correlator  is  to  detereine  output  anplitudes  under 
various  conditions.  The  test  system  is  sho%m  in  Figure  69. 


Figure  69.  Correlator  Output  Test  System 

Simulations  will  determine  output  amplitudes  for  a 
synchronized  correlator  and  for  offsets  between  the  local  and 
incoming  codes.  The  output  amplitude  for  the  synchronized 
correlator  varies  within  a  range  and  will  not  always  match 
the  expected  maximum  amplitude  given  in  Eg  (20) .  This  is 
attributed  to  the  method  of  integration  used  by  the 
simulation  model  and  the  phase  offsets  that  occur  between  the 
local  and  incoming  carrier  frequencies.  The  simulation  model 
for  the  integrators  use  a  trapezoidal  approximation.  The 
results  of  the  simulations  are  given  in  Table  15.  While 
there  is  an  overlap  in  the  output  ranges  between  the  offsets. 


a  clear  distinction  can  be  dram  for  synchronization  within  a 
one>half  chip  period. 

Table  15.  Correlator  Output  Aaplitudes 


Code  Offset  Output  Aeplitude 


(#  Chips) 


Range 


4.20e6 

to 

8.38e6 

0.35e6 

to 

3.80e6 

1.90e4 

to 

1.44e6 

1.60e3 

to 

1.60e4 

1 . 00e2 

to 

2.60e3 
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